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ckapt&A I 


HyPERFIWE WTERACnONS 


A. INTRODUCTION 

Hyperfine 'interactions in general deal with the 
interaction of the nucleus with the charge and current 
distribution of the surrounding electrons and ions. The 
electric charge and current distributions inside the 
nucleus can be expressed by means of a multipole expansion 
of electric and magnetic multipole moments. The interaction 

can be expressed as 

1. The electric monopole term which is the interaction of 

■ * 

nuclear charge distribution with the surrounding 
electrons giving rise to the isomer shift, 

2. The interaction of the magnetic dipole moment with the 
effective magnetic hyperfine field at the site of the 
n ucleus . 


1 



2 


3. The electric quadrupole tensor interacting with the 
electric field gradients. 

The isomer shift results because of the overlap of 
the nuclear charge distribution and the electron density 
at the nucleus. If one makes a reasonable theoretical ' 
estimate of one of the quantities the other can be 
obtained. The relative changes in the nuclear charge 
distribution for different excited states of the same 
nucleus can be measured rather accurately. More often 
one measures the changes in the electron densities at the 
nucleus for different chemical compounds. 

The magnetic dipole moment measurements of the 
nuclear states are possible using external magnetic fields. 
If we know the magnetic dipole moment of the nuclear state, 
one can obtain the information about the induced hyperfine 
fields at the nuclear site due to the surrounding electrons 
and ions. Even the diamagnetic ions embedded in ferromag- 
netic hosts can experience large hyperfine fields. These 
large induced fields arise mainly from the s-type conduction 
electron polarization, core polarization of s-electrons by 
the polarized d or f electrons resulting in a net s-type 
contact spin density at the nucleus, large contributions 
from the unquenched orbital angular momentum and from the 
overlap of the 4s and 5s electrons with the polarized 3d, nf 
electrons of the neighbouring ion. 



Lastly, the electric quadrupole splittings which 
arise because of the interaction of th^e electric quadrupole 
moment tensor with the electric field gradients can be used 
to obtain information about the electronic and ionic contri- 
butions to the electric field gradients at the nucleus. If 
the electric field gradients are known, the electric 
quadrupole moment of the nucleus can be measured. 

Mossbausr effect and perturbed (Y“Y)ang ul ar correla- 
tion (PAC) techniques are used to measure the isomer shifts, 
magnetic splittings and quadrupole interactions. The high 
sensitivity of these techniques make them particularly 
suitable for the studies of the dilute impurity problems 
undertaken in the present thesis. The local environment 
of the probe nucleus plays an important role in all these 
interactions with the result that a proper source prepara- 
tion is very important. In all the PAC measurements 
reported in this thesis, only time -d i<f fere nti al measure- 
ments are carried out and a Fourier transform of the data 
were done so that the observed linewidths m.ay be used to 

infer about the existence of non-unique impurity locations 

0 

in the lattice and the satellite frequencies arising due 
to different lattice locations of the impurity ions. 

It is my earnest hope that the new measurements 
reported here and as well as some of the measurements 
repeated along with the interpretation and discussions 
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given by me extends some of our understanding and/or 
remove some of the existing uncertainties and anomaliss 
about the different mechanisms involved for the large 
induced internal magnetic hyperfine fields in particular, 
and the hyperfine interactions in general. 

B. HYPERFINE INTERACTIONS 

Hyperfine interactions [ 1-4] provide a single area 
of overlap between two rather diverse disciplines of physics, 
viz., the electronic structure of matter and the structure 
of nuclei. These interactions refer to deviations from 
Coulombic laws of interaction between a nucleus and sur- 
rounding electrons. Since a point charge yields, the 
Coulombic force field, the hyperfine interactions thus 
refer to and measure deviations of a nucleus from a 
point charge character. Such deviations are expressed in 
terms of a series of multipole moments. Parity considera- 
tions, then show that the systems with an axial symmetry 
have an alternating sequence of non-zero multipole 
moments, viz. electric monopole (ED), magnetic dipole (Ml) 
electric quadrupole (E2) etc. 

Any hyperfine interaction measurement yields a product 
of two quantities - one characteristic of a nuclear property 
and other signifying the electronic or solid state aspect. 
Thus, an a priori knowledge of one of these parameters is 
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essential as well as helps, determine the other. The 
interaction Hamiltonian of a nucleus embedded in a solid 
can be represented as: 

H = + . . . +] [ + . . .] ( 1 ) 

where represents int ra-n ucle ar interactions, is 

1 Q 

approximately 10 times stronger compared to the remaining 
contributions and is responsible in giving rise to the 
nuclear levels. The other terms, which act as a perturba- 
tion, arise due to the hyperfine interactions. The second 
term in this expansion is the electric monopole term 
which represents the coulombic interaction of the nuclear 
charge and the surrounding electronic charge density. 
is spherically symmetric and has no effect on the m- 
degeneracy of the eigenstates of H|\j . The only effect 
has is to produce a small shift in the energy eigen-values 
of Hj^. Hq represents the interaction of the nuclear quadru- 
pole moment tensor with the electric field gradient at the 
nuclear site. Finally, represents the interaction of 

nuclear magnetic dipole moment with the effective internal 
field at the nuclear site. Both and Hq lift the 
m-degeneracy of the eigen states of . The general level 

s-plit tings are schematically shown in Fig. (1) for the ground 

57 

and first excited states of Fe. 
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1 . I some r Shift 


The net change (Ae) in the Y-ray energy due to 
electric monopole interaction can be expressed as, 

2 


Ae 

!l 

> 

m 

CO 

- A Eg =”ir Ze'^ (0) j 

[<R^> - <R„ 

c g 

>1 

(2) 

w he re 

A E and 
e 

AE are the shifts in 
9 

the energy 

of 

the 


excited and the ground states of the nuclei and and R 

e g 

are the nuclear radii in the excited and the ground states 

2 2 . 

respectively. Since the factor (< R > - <R >) is constant 

B g 

for a given Mossbauer transition, an energy shift would 

occur, only if the source and the absorber nuclei experience 

t 

different electronic charge distribution, i.e. different 
I’l' (0) f . This energy shift is called isomer shift (6) 
and is given by, 

2"^=' ( < Rf>-<R2>,t|^„(D)|2 . 

(3) 


« = = 


e g - - ■ ^ A 

i2 




where the subscripts A and S refer to the absorber and the 

source respectively. If R > R , a positive isomer shift 

B 9 

implies an increase in the s-electron density in going from 

source to absorber. However, if R < R , which is the case 

e g’ 

57 

of -F e , the same p os i t ive is ome r shift shall signify a 
decrease in the s-electron density at the nuclear site in 
going from the source to the absorber. b'e shall defer a 
further discussion on isomer shifts to Chapter IV . 
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2. Electric Quadrupole Interactions 

As the parity cohSiderations exclude the possibility 
of nuclei having an electric dipole moment, the second non- 
vanishing term of the electrostatic part of the interaction, 
is a quadrupole interaction of the nuclear quadrupole moment 
with the electric field gradients produced at the nuclear 
site by the surrounding electronic charges in the crystal. 
The electric field gradient at the nucleus is defined as, 


2 

V V 


3 X i 3 X . 

u 


With a proper choice of axes, this can be reduced to a 

diagonal form, with three principal components, i.e. 

2 2 2 

■ > — ^ ^ . These three components are however not 

3x^ 3y^ 3 2^ 

independent, being connected through Laplaces equation, i.e. 

V^v +’-4 ^0 (5) 

X o y o z 

in the region where the charge density vanishes. It is 
important to note that although the s-electrons do have a 
finite probability of being at the nuclear site, however, 
being spherically symmetric, they do not participate in any 
quadrupole coupling . The electric field gradient is usually 


expressed in terms of asymmetry parameter defined 
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The interaction Hamiltonian is given by. 


H 


Q 



_ e^oQ 
- 41(21-1 ) 



1(1+1 ) + ^ ] 


( 7 ) 

( 8 ) 


The first order perturbation theory yields the energy 
eigen values as, 


E 


m 


I 


e^ Cl q (1 


+ 3 ^ 


f'4 


31 ^ - 1 ( 1 + 1 ) 
t 4I(2I-T1 


(9) 


In ^"^Fe, since the ground state is a spin 72 state, 
it is not affected by quadrupole interaction. The excited 
state however with a spin 3/2 gets split into two levels 
corresponding to nij = +3/2 and mj = +72. 


The quadrupole splitting can be computed using 
equation (9) and accordingly, 


- ^+ 3/2 ~ ^+72 


1 

2 


e^qQt 1 


+ 



( 10 ) 


Thus, we see that a pure quadrupole interaction is 
only partially successful in lifting the m-degen eracy , it 
being a quadratic operator in I^. 


3. Magnetic Hyperfine Interaction 


The first magnetic term of the expansion is the well 
known nuclear Zeeman effect term. The interaction Hamiltonian 
is. 





eff 


( 11 ) 
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where is the magnetic dipole moment of the nuclear 
state and is ihe effective magnetic field at the 

nuclear site. The magnetic perturbation in general 
completely lifts the m-degeneracy of the state. The 
energy separation of the two successive levels is. 


AE^ = < Im[ - y llm> -< Im-1 | - Im-1 > 


“ y ^ I I nij 


( 12 ) 


Since the ground state and the excited state magnetic 
moments are different, the splittings are different. 
Further, the transitions are governed by the selection 
rules. Am = 0,+ 1 (for Fe, e.g.), depending on the 
nature of the transition. These two factors coupled 
together are responsible in giving rise to a six finger 

57 

pattern in Fe. The effective magnetic field seen 

by an impurity in a ferromagnetic host like Fe, Co and Ni 
was first discussed by Marshall [5]. He gave 


»Bff ' "ext - “Ns + X «s 


H 


hf 


(13) 


H , , DM and — M are the external, demagnetizing and 
the Lorentz fields respectively. The origin of these 
contributions are well understood and can therefore be 
well accounted for. These are typically of the order of 
a few kOe. Demagnetization factors for various source 
shapes have been tabulated by Bozorth [6]. The saturation 
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magnetization is 7, 6, 2 for Fe, Co and IMi [7] respec- 
tively. H is a small correction to the Lorentz field 
for materials having a symmetry other than cubic* Even 
for hexagonal symmetry (e.g., Co at room temperature), 

“j" ^ 

turns out to be about a thousandth of the Lorentz field 
[8] and hence can be altogether neglected. The last term 
represents the true hyperfine field and we shall discuss 
this in detail. The various contributions to can be 

given as [91, 


'^hf = "ext + 


orb 


+ H , . + H + H ^ + H ■ (14) 

dip core cond vol 


Horb ^spresents the hyperfine field produced by the orbital 

current of the electrons and can be given as, 

-*■ 

l.-I 


H 


orb 


2^ < 

I . ^ B T 
1 r . 


( 15 ) 


. th 


v-jbere 1^ represents the angular momentum of the i- 
electron having a position vector r^. Due to spherical 
symmetry, the orbital contribution is zero for the closed 
shells and for the half filled shells. Also, the orbital 
contribution is small when the orbital angular momenta is 
quenched by the crystal field gradient of the surrounding 

-i- 

charge distribution. ^^iip ® spin dipolar 

interaction between the nuclear magnetic dipole moment and 
the magnetic moment of the electrons. 


- -t^) ^r^^T) - r? (s .1) 

"dip = I < > 

1 r . 


( 16 ) 
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where represents the spin of i electron with r^ as 
its position vector. This again vanishes for system with 
a cubic symmetry. 

The other two contributions are the core polarization 
and the conduction electron polarization arising due to 
spin polarization of the core and conduction band electrons 
respectively. The exchange interaction between the 3d{4f) 
electron spins and the core s-electron spins causes a 
'repulsion' of electrons with antiparallel spin and an 
effective ’attraction’ with the electrons having a parallel 
spin. Such an interaction therefore results in a difference 
in the radial wave functions for electrons with opposite 
spin, thus resulting in spin polarization of the core, 
leaving a net down spin density at the nuclear site. This 
induces a negative hyperfine field at the nuclear site 
via a Fermi contact term given as, 

I 1 I I (17) 

ns - 

A similar exchange polarization of conduction elec- 
tron spins is also produced via an exchange with host and 
as well as with the solute d ( f ) -electrons , which again 
results in a net hyperfine field at the nucleus. 

Lastly, ® positive field arising due to 

valence ns-electrons , when the volume’ of the solute atoms 
is larger than the volume available to it on removal of a 
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host atom. The ns-1 
matrix and get posit 
tional to the voluma 
contribution is howe 
reasonings. Chapter 
and the relative mag 


ike electrons overlap with the host 
ively polarized by an amount propor- 
misfit of the solute atom. This 
ver based somewhat on empirical 
V deals in detail with the origin 
nitudes of these contributions • , 
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chapte.A II 


EXVERTMENTAL TECHNIQUES 


A. INTRODUCTION 

The methods employed in hyperfine interaction studies 
can be classified [1] into two categories depending upon 
whether the impurity is stable or radioactive.. The methods 
used for stable nuclei are; 

(a) Hyperfine structure of optical spectra 

(b) Atomic and molecular beam methods 

(c) Nuclear magnetic resonance 

(d) Electron spin resonance 

(e) Nuclear specific heat method 

(f) Interaction of polarized neutrons with polarized nuclei. 

The methods employing radioactive nuclei are: 

(a) Mossbauer effect 

(b) Perturbed angular correlation and its variants 

(c) Nuclear orientation. 

15 
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Besides these techniques, quite a few more methods 
employing hybrid principles of the above methods have been 
used. A typical example is f\jHR/PAC, the resonance destruc- 
tion of angular correlation. The perturbed angular correla- 
tion and iMossbauer effect techniques used for the measure- 
ments presented in the present thesis are described in 
greater detail in the following sections* 

B. PERTURBED ANGULAR CORRELAT IONS 
T. Introduction 

In I 95 O 5 Brady and Deutsch [-21 recognized the possi- 
bility of measuring the g-factors of excited nuclear states 
using an anisotropic angular correlation between gamma rays 
emitted in cascade. The first such measurement was reported 
by Aeppli et al. [3] . Since then the theory of extra- 
nuclear perturbations on angular correlations have been 
developed to a high degree of completeness and has been 
dealt with extensively in reviews by various authors [4-11]. 

In particular, the applications of PAC to solid state problems 
have been discussed in detail by Heer and Novey [ 9 1. More 
recent reviews by Shirley and Haas [10 ] and by Deutch [11] 
discuss at length the present trends in PAC rose arch , F ig . (1). 
In the following a brief and a self contained account of the 
theoretical formulation, relevant to the measurements presented 
in this thesis is given. 
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Schematic’ ■chairact.Erization of PAC experiwsots b.as®:'^' -oii C-a) 
of praceeding decoy , (b) Physical state of tha^> ample and 
(p) Type of extra' nuclear perturbation, Cach box is described 
three coordinates, corresponding to different permutations of 
the above three factors, » 
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2. Theoretical formalism 

In a radiative nuclear decay, through a three level 

system, via successive emission of two photons, an angular 

correlation in general exists between their (photons) 

propagation directions. Thus in an angular correlation 

experiment, one looks for the probability W ( , k 2 ) d Q ^ d fi 2 ? 

Rl , R2 

that a nucleus decaying through a cascade ^ I ^^f’ 

emits two radiations and R 2 in the direction and k 2 
into solid angles d and d respectively, where W(k^,'^ 2 ) 

is the angular correlation function. For a system in a 

# 

field-free space, i.e. for isolated nuclei, this correlation 
depends upon and yields information regarding the level 
spins and the transition multipolarities of the radiation 
fields. This angular correlation, however, is perturbed, if 
the intermediate state interacts with the environment. 
Physically, an extra nuclear perturbation lifts the 
m-degeneracy of the intermediate level and changes the 
m-sublevel population, consequently altering the angular 
correlation pattern. 

Any angular correlation experiment can be described 
using the following three principal tools 112-143. 

(a) The Density Matrix Formulation, which provides a 
conceptually simple and a mathematically elegant 
method to describe an incoherent mixture of impure 
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or mixed states, which is the case in a nuclear process 
like angular correlation. Impure states arise due to 
the fact that in an angular correlation experiment, 
we are often met with transitions of mixed multipola- 
rities and that not all the quantum numbers are known. 
Incoherence in these states arise due to the fact that 
in a nuclear process, which is essentially a statistical 
phenomenon, each member of the nuclear ensemble develops 
in time in a statistical manner, totally independent 
of others. 

(b) The Rotation Matrices, the usage of which gives one 

a freedom to calculate transition matrix elements in a 
convenient quantization frame of reference and then 
finally rotate them, with the help of rotation matrices, 
so that they correspond to actual experimental conditions, 
lastly 

(c) R ac ah ' s An gul ar Momentum Algebra, which gives a 
convenient vjay to perform summation over various 
m-substates. 

i. The unperturbed correlation function 

A typical angular correlation cascade is shown in 
Fig, {2a) . The initial nuclear state ll^m^>, described 
totally by the density matrix P^, decays by emitting 
R^(L^,l!j) in the direction k^, to the intermediate state 




(a) Schematics of. Unperturbed Angular Correlation 




(b) Angular Correlation 

^ 'Scfteiaaiic's of ( ) ^'inyular Carredationi ar>rl (b) T'erbui: 

Angular C Dxr el -j I: ^ , ' ■ 
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} Im > de scr ibe d by p (k^). The intermediate state in turn 
decays by emitting R2(L2,L2) to the final state ] > 

described by p(K-jsK2). P(k^>k2), thus contains all the 
information about the system and the radiation pattern. 

If and H2> be the interaction Hamiltoniany causing the 

decays and R2, we have. 


< Im I P(k^)l ]m'>= I <lm| I^rn^x I^mJ ^ ±\ ^ > 

m m ^ 


m . , ml 

X X 


< Im’l H, ! I .m'. >* 
1 1 1 


5^ denotes a summation over all the unobserved properties 
like, e.g. the polarization. Similarl-y then for R 2 » we 
can write, 

< I^m^ I p (k^ ,k 2 ) I I^m^> = S 2 I < I^m^l H 2 i Im xlm 1 ftk^)\ 

m, m' 

Im‘>< I^m^‘ I H 2 I Im' > 

Using equation (1) we have, 

< I^m [ P( k^ , k 2 )l I^nij,> = 5^52 I , I < H 2 1 Imxlmj 1 

m . , m . m , rn ’ 

I.m.Xl.m.I P . 1 I .m’>* < I -m'l |lm’> * (2) 

xxxxxxx fr^ 

For an isotropic initial state, the density matrix is 
diagonal , i. e . 


< I .m. I p. |l .m’. > = (21 . + 1) 

11 i ix X 


-1 5 

m . m . 
X X 


The correlation function is given by. 
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W(1<^,k2) = I< P(k^,k2)l > (4) 

mf 

where < [ p(k^ >k 2 ) I > is the probability of finding 

the nucleus in the state I I^ni^> . Thus substituting equations 
(2) and (3) in equation (4) we get, 

Vi/ ( k^ , k 2 ) = 5 ^ 52 I I I < I^rn^l H 2 I Im >< Inj H^( I > 

m- m . m , m ’ 
f 1 

<Ijn'| H.J H 2 I I m'>' (5) 

This is the unperturbed correlation function or the 
true correlation function, i.e. nothing happens to the 
intermediate nuclear spin, during the lifetime t of the state. 


ii. The perturbed correlation function 

An extra-nuclear perturbation, in general, lifts the 
m-degeneracy of the intermediate level and changes the level 
densities during the lifetime T . The density matrix p {1<^ ) , 
therefore changes in time (Fig. 2b). Thus to obtain an 
expression for the perturbed correlation function, the time 
evolution of the density matrix has to be introduced in the 
formalism. The time evolution of the density matrix is 
described by the Heisenberg equation of motion, i.e. 


9p 

3t 


-- [k , p] 

K 


which for static interactions has the solution, 


-1 


p (k^,t) = A“ (t) p(k-j,t=:0) A(t) 



where 1C is the interaction Hamiltonian, describing the 
extra nuclear perturbation and 
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A ( t) 


-iict/K 

e 


The perturbed correlation function is then obtained 
by replacing p ( kp by p(k^,t) .Fig, (2 b)). Thus we have, 

P ( k.j ,k 2 ,t; I > = 5^ I H^I ImJ A(t) 

m , m, 
a b 

p(k.,,t=D) A“ht)l Im|^> <I^m|lH2l 

= I I < I.pm^jH 2 l Im Xlm ] A(t)l Im > 

rn g y m y m 

/^b 

< Im I p t it.] , t=0) 1 Im'XIm’ 1 A“^ ( t) I Im^^^ 

Using equations (1) and (4) and once again assuming the 
density matrix to be diagonal in the initial state, the 
perturbed correlation function can be 'written as, 


W( k.] , k^t ) 


I I 

9 m, m’ 


I < I^m^ H^l Im^Xlm | \ 


< Im' I H] 


* 

I . m •> 

1 1 


< I^m^t 




< Im^ 1 A( t)|lm><Imj^lA(t)llm’> (7) 


A comparison of equations (5) and (7) indicates that the 
effect of extra nuclear perturbation is contained in the 
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matrix elements, 


Im m l G(t)l Im m'> = '^Im MMt) { [A (t) | Im’>*{8) 

a b a b 

Using rotation matrices and the Racah algebra, the 

P AC function can be written in the conventional form, i.e. 


W(k^,k2,t) = ^ 


1 5 ^2 

InI ^ y l\l2 


/TTk^jTTTTTi^^ 


N N N ^ N. ^ 

Sk,i<2 <h' 

whare the ^rotated* perturbation factor is related to the 
original perturbation factor through the relation, 


^2 


2l4-m +m, 

I (-) “ ^ 

m , m 


[ (2k^-Hl ) (2k2+1 )] 


Ilk. I I k^ 

) ( ^ 

m* -m N . m, - m 

1 b a 2 


) < E(t) jlmm*> (10) 


In the above L and L* are the angular momenta carried 
by the y-rays. For mixed multipole transition , both L and L* 
are non-zero. The summation indices k and N are determined 
by the selection rules IL-L*! ^ k L -h L [ and iM^ = m’ -m and 
Np = m. -m_ and the coefficients, 


^k^ ’*”1 ’ h’ 


1 ‘-1+^i 

-i-2 (L^,L^ ,L.,I) + (-) 26 


F^(L^,L;,I.,I) + 6^ F^ (L',Lpi.,I) 
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and the mixing ratio 

<I {[ L-l ^^11 I > 
<I II TT^I! I > 


s = 


The F-c oeff icients are tabulated by Frauenfelder and 
Steffen [4] . Chosing the axis of quantization to be the 
symmetry axis, we have. 




y2 I I k 

(t) = I [ (2k +1)(2k +1)] ( ‘) 

m m’ m' - m N 


( 


I I k 

m ' - m l\l 


^ -i(E -E„, )t/)( 
2 . mm' ' 

) e 


(11 ) 


For static magnetic perturbation, 

Ej^-E , = < Im! - J j.’h I Im> -< Im-l| -•Uj.H[lm-1>= K W|_ 




the perturbation factor reduces to a simple form^ 

y2 -iW| t Ilk Ilk 


( 2 k^+ 1 ){ 2 k 2 + 1 ) ] e 


I ^-m m’ N ^^-m m' N ^ 


m 


-iw t 
= e ^5 


ki t<2 


( 12 ) 


Substituting (12) in ( 9 ) we get, 

-*■ ^ -iw, t -> 

W(k,, ,k 2 ,H,(t) = I (R.,) (R 2 ) e ^ Yj^ ( k^ ) 

k-| k2 1 1 1 

^k2 ^*^2^ 


(13) 


■ The ’observed* correlation pattern is finally obtained 
by multiplying the above equation by the probability that the 
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radiation will be emitted at all, at time t, after the 

emission of R-j . This is given by a Poisson distribution 
1 -t/x 

— e , where T is the mean life of the intermediate state. 

T 


Equation (13) describes the correlation pattern that 

would have been observed experimentally, if the correlation 

apparatus had a ’perfect' time response. In practice, however, 

as the experimental system has a finite resolving time, the 

observed angular correlation pattern depends on the relative 

magnitudes of the resolving time x ^ of the detectors and 

the lifetime of the intermediate state x . Thus, if X »x, 

o 

one can then measure only a time integrated correlation 

pattern (IPAC), i.e. 

00 


W( s,<», H) 


i W(e, t, H) 


dt 


“ -■t/x,, 

e d t 


(14) 


which is independent of the apparative constants. On the 
other handj if the apparatus has a resolving time such 
that T ^ << we get a time differential pattern (TDPAC), 
io3. the pattern can be described by equation (13) and 
hence one is able to observe the rotation of the angular 
correlation pattern with time. Thus ^ a TDPAC measurement 
directly 1 senses’ the motipn of the nuclear spin as a 
function of time ^ i,e. 


W(6 , t,H) = 




dt 


'1 




where min ( t_-t. )>> x . 

2 1 * o 
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3. Experimental geometries 

It is often possible to select the propagation 
direction in such a way that facilitates the extraction 
of the perturbation information from the experimental data. 

There are three variants of PAC set-ups employed mostly. 

i. Perpendicular geometry 

In this most widely used geometry, the magnetic 
field is perpendicular to the plane of the detectors. Thus, 
referring to Fig. (3), 

0^ ^ = tt /2 and 0= <j) ^ 2 

In such a case, the correlation pattern reduces to 

a particularly simple form, 

k 

max 

W, ( , t, H) = I Aj^ (R 2 ) Pj^Ccos (e + W|_t)3 (15) 

l <=0 

even 

or equivalently, 

k 

max 

Vi (9 , t , H ) = I b| c os k { ® ) (15) 

-L k=0 ^ 

even 

Interpreted on the basis of vector model, equation 
(15) offers a simple picture, that is, that a magnetic 
interaction results in a Larmor precession of the nuclear 
spin by an angle w^^t and thus al-ters the angular correlation 
pattern appropriately. The direction of precession is determined 
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by the requirement that the energy E = and the torque 

T = U X H . Figure (4) elucidates this point clearly. Thus, 
in a fixed angle (9) geometry, the experimental data is a 
decay curve superimposed with a harmonic part. The modula- 
tion frequency in such a time differential method is 2wj^. 

The method further allows a determination of the sign of the 
field, by observing the sense of rotation. 

The time integrated correlation (IF AC) in such a 
geometry is given by 



It is therefore obvious that in an integral perturbed angular 
correlation, the angular correlation pattern is both rotated 
and attenuated. Further, it is worthwhile to note that a 
time differential measurement doss not require an a priori 
knowledge of nuclear decay parameters, while an IPAC measure- 
ment does assume an a priori knowledge of these parameters. 
The perpendicular geometry DPAC was used in the measurement 

of hyperfine field on 5cC o and'Taiy: and in the g-f acto r 

181- 

Ta. 



measurements on 
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ii. Randomly oriented domains technique 

It was pointed out by Matthias et al=-El5] , that P AC 
could also be observed in a randomly oriented domain sample, 
i.e. when no external polarizing fields are employed. Here, 
therefore, although, 

8 ^ + 0 2 = TT snd 0 = <i) ^ _ 4>2 = 

are fixed, however, 6^ may take any value. Accordingly 
then to get the correlation function, an integration over 
all the possible values of9j_, has to be carried over. The 
correlation pattern in such a case is a superposition of 
k frequencies, weighted equally, i.e. 

w (e. H, t) = I A^R,) I cos' NW|_t 

k N 

P|^(cose ) ] (17) 

which for k=2, takes the form, 

VJ( , H, t) a [1 +2cosw^^t + 2cos2w^t] 

The integral correlation pattern, however, remains 
the same as in perpendicular geometry. This TOP AC/randomly 
oriented domains method was used in hyperfine field measure- 
ment on Ta in Co matrix. 

iii. Raghavanand Raghavan’s geometry 

Raghavan and Raghavan E 16^ pointed out two geometries 
that emphasize a modulation frequency of W|^ rather than 2wj^. 
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The method, however, is not capable of giving the sign of 

the f ic Id „ This geometry is characterized by, 

# 

The polarizing field now lios in the plane of the 
detectors. This technique was used in the measurements 
on Re in Ni (not presented in the thesis) . 


4. Data Reduction 


In the data reduction the first part is to extract 
the frequency content of the observed TDPAC spectra. This 
is easier after correcting for the exponential decay of the 
nuclear excited state. This in most cases is accomplished 
by forming a reduced counting rate ratio R(t) defined as, 
e.g. , for the perpendicular geometry one has, 


W(e,t . , H+ ) - W(e , t . , Hi) 
R (6 , t . H ) = 2 ^ ^ 


3 sin 2 w, t . 
22 L 1 


W( 6 ,t^, H+ ) + W(e , tj_, HT) 


4 + A 


22 

(18) 


The two correlation patterns corresponding to the 
Hf(up field) and 11 + (down field) are phase shifted with 
respect to each other by tt. 

Equivalently in Raghavan and Raghavan's geometry one 
calculates the reduced counting rate by computing the ratio 
W(ir) - W(it/2) 


R( t) 


I A22 cos w^t 


W('ir) + W('rr/2) '’l + 1/16 A 29 + 3/16 A 22 cos 2w^t 


for A 22 « A 44 


(19) 
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Alternatively, one could also form a ratio C(t) defined as, 


C ( t ) 


W 




A 


H. t) - 

-t/ T 

e 


A 


-t/ T 


which would eliminate the exponential decay factor. 

From these ratios giving the reduced counting rates 

R(t) and C(t)'s, the frequency content could be directly 

computed by an extremely elegant and powerful method, 

suggested by Matthias et al. [17], who observed that 'except 

for the details of the line shape, the fourier transform 

of the correlation spectra, decay part taken out is 

essentially analogous to the free precession version of an 

— 1 ' 

1\1MR line (i.e., no R.F.) with a line width T . Thus to 

extract the frequency content, one just has to calculate 

the Fourier transform of the, reduced counting rate defined as. 


00 


A ( w ) 

= / R(t) 

sin 2tt Wj^t 

dt 


o 



00 



B(w) 

= / R(t) 

cos 2tt w^t 

d t 


O 


F(w) 

= [ (A^(w) 

+ B ( w 



( 20 ) 


Lastly, the range of the applicability of a TDP AC 
experiment is determined by the relation 




T <. T 

Wj^ 


represents the time period corresponding to the 
w^ 


where x 
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Larmor frequency. The presently available time resolutions 
give the lower limit of 1 nsec and accordingly TDPAC is 
applicable to the lifetime range 


-9 -2 

10 s < '.T <10 s 

The upper limit however for conventional off-beam 
coincidence methods is 1 sec, because of poor true to 
chance ratios. However, the upper limit can be considerably 
larger for in-beam experiments. The uncertainties in most 
favourable cases are 


C. M0S5BAUER EFFECT 
1 . I ntro duction 

The discovery of recoilless gamma resonance fluore- 
scence by Mossbauer [ 1 8 ] brought with it a diversity of 
applications, diffused practically in all the domain of 
natural sciences. Although, in principle, the effect 
originally belonged to nuclear physics, its solid state 
applications are now dominant, to the extent that it has 
become one of the standard tools in the investigation of 
solids. Such a multitude of applications have arisen 
essentially due to extreme energy precision offered by 
Mossbauer effect with relative simple experimental equipment. 
Probably no other technique offers so much of a v^/ealth t)f 
information with so great an ease. Various aspects of 
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Mossbauer effect, its principles and applications, have 
been extensively reviewed [ 19-24] and accordingly we shall 
briefly outline the underlying principles and indicate the 
variety of information that can be gained by a typical 
Mossbauer study. 


2. Theory 

Mossbauer showed that the nuclei that are embedded 

in a solid can emit and absorb low energy gamma rays in a 

recoil-free manner. A free nucleus, while deexciting from 

an excited state E^, by emitting a gamma ray of energy E_^, 

2 2 

recoils with an energy , where m is the mass of 

the recoiling nucleus. Conservation of energy gives the 

gamma ray energy E^ = E^ - Ep. Analogously for the absorption 

process, it follows that in order to excite nucleus from 

its ground state to the same excited state E^, a gamma ray 

with an energy E = E + En is necessary. Thus, there is 

Y o K 

a net deficit of 2 Eq, for the resonant absorption of the 

gamma ray. Resonance absorption can still occur if 2Ep<. T, 

the width of the nuclear state E (determined by its life- 

o 

tims). Typically for ^"^Fe, E^=14.4 keV, 2x10 ^ eV 

-9 

and r~ 5 X 10 eV and obviously enough it is impossible to 
get a resonance under such conditions. Things, however, are 
different when the nucleus is rigidly bound in a lattice. 

The recoil momentum is now taken by the' lattice as a whole . 

In such a case, the .inevit ab le reco il energy is drastically 
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2 ^ 

reduced, and is now given by Ep= Ep/2 Me" v^/here M is the 
mass of the lattice instead of m, the mass of the nucleus. 
Under such conditions the recoil energy is not adequate 
to excite the solid to. any higher phonon states so that the 
condition 2E^<<r is amply satisfied and accordingly resonant 
emission and absorption occurs. Here, in principle, most 
of the emission and absorption processes occur recoil-free, 
i.e. the fraction (f) of the gamma ray emitted recoil free 
can be almost close to one, while for a free nucleus f =0 . 

A real crystal however does not provide an infinite binding 
to a nucleus and accordingly the internal degrees of freedom 
of the lattice may be excited by the emission and the 
absorption processes. In these cases, f may be between the 
two extremes, i.e. 0^ f <. 1 , depending on the energy of 
the transition and the 'rigidity' of the lattice. This 
may be summarized by the Lipkin’s sum rule [25] 


I 


[ E„ - 1 P(n.-n ) 

Of n. f X 


(21 )■ 


where E and E are the energies of the lattice in the 
n . n „ 

1 f 

initial and final states ( n^ S. n^) and P(n^,n^) is the 
corresponding transition probability, representing the 
fraction of events inducing a lattice transition from a 
state n^ to state n_p. Clearly then, the transitions of 
interest in Mossbauer effect are identified by the 
probability P ( n^, n^), which is, in fact, the recoil free 
fraction. ForaDebye solid [ 4 3 , 
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P(n.,n.) 


= f 


2 2 

exp [ -k <x » ] = exp [ ■ 


2 


TT 

8g2 


) ] 


( 22 ) 


where 0g is the Debye temperature and k^ is Boltzmann 
constant. Thus, the conditions necessary to observe recoil 
free transition and hence the Mossbauer effect are given by 
equation (22); 1) and i. e . , should be small, 2) lattice 

temperature should be small, 3) the Debye temperature of 
the solid should be large. 


To summarize, therefore, under these conditions ,. the 
emitted gamma ray has the full energy E^ of the nuclear 
transition and has a natural width (r ) (which is indeed the 
key to high energy resolution offered by Mossbauer effect) 
and accordingly have the right energy to be reasonantly 
absorbed by a similar nucleus under similar conditions. One 
is then able to study small perturbations in the nuclear 
energy levels arising due to, e.g., hyperfine interactions, 
by observing the change required in the gamma ray energy 
for it to be resonantly absorbed. Such changes are 
afforded by an appropriate Doppler shift. 


A typical Mossbauer spectrum is characteristic of 
the atom under study and depends critically on the valency 
and the symmetries of the surroundings of Mossbauer atom. 
When several phases are present, the resultant Mossbauer 
spectrum is a linear superposition of the constituent 
spectra of differbnt phases. Further, for thin samples. 
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the intensity of Mossbauer spectrum is proportional to 
the concentration of Mossbauer nuclei. In the followingj, 
we shall briefly discuss the information that could be 
derived from a Mossbauer spectrum. 

A study of the nuclear hyperfine splittings directly 
reveal nuclear electromagnetic coupling coefficients both 
in the ground state and in the excited state. Then, the 
isomer shift, which can be measured only by Mossbauer 
effect is given by the shift in the centre of gravity of 
the spectra with respect to zero velocity. The isomer 
shift directly reflects on the s-electron density at the 
nucleus and consequently yields information regarding the 
chemical bondings and site symmetries. There is yet 
another contribution to the centroid shift, which is 
strongly temperature dependent, viz., the second order 
Doppler shift. The second order Doppler shift could be 
understood as follows. The relativistic energy of the 
source as seen by a moving absorber is 



1 ■ + 


;v 

c 




arising due to the fact that both the emitting and absorbing 

nuclei are vibrating inside the lattice with very high 

12 

frequencies ("10 Hz). In such cases, although the 

average displacement during a spectroscopic event {— term) 

2 ^ 

averages out to zero, the term — however, does not. This 

c 
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results in a small perturbative addition to E and accord- 

o 

ingly results in a centroid shift, depending upon the mean 
square velocity of the nucleus. The temperature dependence 
of such ,a term is obvious. 

The next information which once again relates to 
lattice vibration is the recoil free fraction. It was 
seen in< equation (22) , that for a Debye solid, where all 
the nuclei occupy equal sites, the recoil free fraction 
f is 

f = exp (-l<^< x^> ) 

where k is the wave vector of the nuclear radiation. The 
recoil free fraction of a nucleus is directly related to 
the area of a Mossbauer spectrum corresponding to the 
resonance absorption, i.e, 

A = f / [ 1 - ]dv = f / [1 ~ dE 

—00 ^ * —oo 

where (1 - ) simply denotes the intensity of absorption. 

f and f' are tho recoil-free fraction in the source and 

the absorber. n is the number of nuclei/cm'^ in the absorber. 

The inference about the recoil free fraction is drawn 
from the temperature dependence of the area of the spectrum. 
It is important to note that the recoil free fraction is. 
dependent upon the mean square displacement of the nuclei, 
while the second order Doppler shift is a measure of the 
mean square velocity. Indeed therefore the f-f actor and 



40 


the second order Doppler shifts provide complementary ways 
of studying the lattice vibrations. The former is dominated 
by the low frequency vibrations whilst the latter depends 
on higher frequency modes. 

The above expression for f is valid only if the 
lattice vibration of the Mossbauer nuclei are isotropic. 

In such cases, f is independent of the propogation direction 
of the radiation. An anisotropy of vibration can exist, 
which in turn results in an anisotropic recoilless fraction. 
For an axially symmetric case, it is given by [21] , 

f (6 ) = exp [ -k^ < X^> ^ ^ ^ p ^ ^ ^ (24) 

where 9 is the angle between the symmetry directions and 
the direction of propogation of the gamma rays. Such an 
anisotropic f leads to an asymmetry in the intensity of 
quadrupole splittings in polycrystalline samples. 

Conversely therefore an unequal intensity of quadrupole 
lines may reflect the anisotropy of f and hence about the 
crystal binding, which is known as Goldanskii - Karayagin 
effect . 

Intensity differences in different hyperfine compo- 
nents can arise from two other sources as well. The first 
one is due to the difference in the relative transition 
probability of different hyperfine components and the 
second one is due to the difference in the angular dependence 
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of radiation intensity of different components. This, 
however, averages out for polycrystalline samples. Thus 
the intensity difference of different hyperfine components 
useful information regarding the relative transi- 
tion probabilities. 

Lastly, the information regarding the electronic 
spin relaxation effects can be inferred from the line 
shapes and line widths. 

Thus, it is by now very clear that every bit, every 
minutest detail of a ["lossbauer spectra carries significant 
information and this indicates that an extremely careful 
analysis is absolutely essential for exploiting the latent 
information. 

3. Experimental 

In Mossbauer spectroscopy a precision Doppler velocity 
spectrometer is required to produce small shifts in the 
gamma ray energies. The nuclear absorption at resonance 
is detected in terms of decreased transmission of gamma • 
rays through the sample containing the absorber nuclei. 

Now, since the natural line width of the 14.4 keV 
57 _g 

excited states of Fe - 5 x 10 eV , corresponding to a 
velocity- 0.1 mm/sec, the accuracy of Doppler velocity 
required is much higher, so as to avoid any instrumental 
broadening of the spectral lines. This puts stringent 
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requirements on the frequency and the amplitude stability 
of the Doppler drive signal. In the following, we shall 
describe the sinusoidal drive spec t r cmeter [2 6] which was 
used for the measurements presented in this thesis 
(Chapter IV ) . 

Sinusoidal drive spectrometers are probably closest 
to an ideal system for two reasons [ 27 3, ( i) that they do 

not put any mechanical stress on the transducer, ii) with 
a sinusoidal waveform, one could measure even small isomer 
shifts giore accurately compared to the constant acceleration 
systems. The only disadvantage of such a drive system is 
that the spectra obtained have to be linearized since in 
this case, the time spe nt/c hannel is not constant. However, 
such a linearization can be afforded by a simple computer 
programme. 

A schematic block diagram of the spectrometer is 
given in Fig. (5). The drive signal was derived from an 
100 kHz, crystal controlled oscillator so as to get an 
frequency stability of one part in 10^. The 100 kHz 
signal after proper level shifting is scaled down by 10 
to get 10 kHz output which is used to advance the memory. 

For such channel advance pulse, the velocity reference signal is 
required to have a frequency (10/512) kHz for a 512 channel 
MCA. This was achieved by scaling down the 10 kHz output 
by 512 ( = 4 X 32 x 2 x 2) using TTL-JK flip-flops. The 
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last stage in this chain is deliberately taken to be — 2 , so 
as to get square wave at the output of the divider chain. 

The reference sine 'wave for the velocity drive signal is 
finally obtained by extracting the fundamental component 
of the square wave through an active filter. A square 
wave input to the active filter minimizes the distortion 
in the sine wave as the square wave does not contain any 
second harmonic. In this case the rising and falling edge 
of the square wave then corresponds to zero of the sine 
wave and hence to zero velocity. Thus if multiscaling is 
started using this sine wave the initial channel {zeroth 
channel) would automatically correspond to zero velocity. 

It is however preferable to have the display with zero 
velocity in the centre. This in turn requires a delay of 
T/4 between MCS start and the filter input. Such a phase 
shift of /2 is achieved by triggering a flip-flop using 
a complementary output of -s-Z counter. 

A high amplitude stability in the amplitude of the 
sine wave is achieved by using a variable gain amplifier 
with an automatic gain control. The output of the amplifier 
is precision peak detected and is fed to the error amplifier. 
The error amplifier compares this voltage, with the reference 
voltage and the difference provides the AGC voltage. When 
the peak of the sine wave goes above or below this reference 
d-c voltage, the gain of the amplifier, is appropriately 
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changed. The amplifier is a difference amplifier, whose 
gain is ccntrclled by varying the base voltage of the 
current source. 1 his peak stabilized - sine wave was 
used as a velocity reference signal. This signal is fed 
to the drive coil of the transducer through a power ampli- 
fier. The transducer converts the electrical signals into 
mechanical motion via the drive coil. A second coil, the 
pick up coil (we used Sanborn 6L\/1 coil) then measures the 
instantaneous velocity of the system. The pick up signal 
is generated by the motion of a pencil magnet rigidly 
connected to the drive coil and this moves in the centre 
of the pick up coil. The pick up signal is compared with 
the reference signal by an error amplifier and such a feed 
back configuration ensures that the actual velocity waveform 
of the transducer is the same as that of the input signal. 
The error amplifier used was a push-pull type class AB 
amplifier using a complementary symmetry circuit. The 
above circuit was found to have an amplitude stability of 
better than 

4. Data Reduction 

Since we scan the full sine wave in 512 channels, 
we get in the MCA memory, two spectra, mirror images of each 
other. The symmetry lies around the folding point, which is 
a velocity extremum point and canno be c alculate d . in a 
straightforward manner. Now, as the full wave corresponds 
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to 512 channels, the zero velocity point would be 128 
channels off from velocity extremum point and hence an 
addition and/or subtraction of 128 channels from the 
channel number corresponding to the velocity extremum 
point would yield zero velocity point. In this way, the 
zero velocity is unambiguously fixed up without using 
any external calibration standard. The zero velocity 
fixed in this way compared favourably with the zero 
velocity obtained using standard sources and absorbers. 

The channel number counted from the zero velocity channel, 
is then linearized by a simple programme. 

D. SOURCE PREPARATION 

The basic requirement in dilute impurity hyperfine 
interaction studies is a dilute source with impurity mostly 
in substitutional sites. . This indeed requires 'a careful 
source preparation. 

One of the standard and fairly reliable methods of 
source preparation is by alloying. The usual criterion 
is that the solid solubility of the solute in the host 
should exceed the impurity concentration desired. The 
solubility is usually good, if the ionic volumes are not 
too different. A rough-guide to the solid solubility is 
provided by the Darken-Gurry plot (Fig. 6). However, 
alloying is not feasible when the melting points of the 
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salute and the host are too different. The other alter- 
native method of source preparation is by high temperature 
diffusion, after electroplating the impurity on the host 
metal. Yet another method is a brute force technique when 
grinded samples are thoroughly mixed up and then solidified 
under a high pressure with a subsequent annealing. In 
such methods, however, the impurity position is always 
quite uncertain and a variety of lattice sites may be attained. 

One of the recent and most promising techniques of 
source preparation is by implanting the. impurity atoms 
into the host matrix by means of an isotope se ap arat or- ' [ 2 8 ] 
or through a recoil after a nucl e ar ■ re action or Coulomb 
excitation. Here, one can, in principle, and under 
suitable experimental conditions can implant any impurity 
in any desired lattice. Subsequent ion channelling 
experiment can help to determine the exact impurity location 
in the host lattice. 
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PILUTE IMPURiry f/yPERFIWE EIELP 
MEASUREMENTS USING PERTURBEV 
ANGULAR CORRELATION TECHNIQUE 


A. HYPERFINE FIELD MEASUREMENTS ON SCANDIUM IN COBALT MATRIX* 
1 . Intraductio.n 

The spin rotation of 1"'', 68 keV state was used to 

measure the hyperfine field on scandium as a dilute impurity 

in cobalt matrix. The 78 kBV-68 keV ( Y~Y) cascade populated 

44 44 

i-n the decay of Ti to ' Sc was used for this measurement. 

44 . . 

The relevant part of the level structure of Sc is given in 
Fig. (1) [1]. Using time differential perturbed {Y“Y) 

angular correlation technique, we earlier reported [2] the 
g — factor of the 68 keV state to be* g = 0.342 +_ 0,002. The 
half-life of this state was measured to be Ty^^ 158.8+2.0 
nsec [2 ] . 

2. Source Preparation 

For the hyperfine field measurement a TiCo alloy 
*P ubl ishe d in P hys . Lett. 47 A , 31 9 { 1 9 74) . 

l.l.T. i . r-FUR 
( ^f NTRA. L UBRARY 
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containing tracBS of Ti activity was preparad. The ^^Ti 
activity, produced through a 5c{p,2n|'^Ti reaction, was 
procured in a solution form (as chlorids) from Nuclear 
Science and Engineering Corporation (USA). The liquid 
activity was evaporated on to a thin cobalt foil ( 99 . 99 %). 
This foil was then melted in an inert argon atmosphere, 
using a radio frequency induction furnace. The molt was 
then gradually cooled down to room temperature. The 
resulting shining ball was rolled into a thin foil, which 
was then annealed for 24 h, and finally cooled down to 
room temperature over a period of 3 h. 

The solid solubility of Ti in Co has been reported 
to be about~ 8 at ^ [3], The atomic radii of Ti and Co 
are 1.47 % and 1.25 1 \ respectively and therefore we expect 
that we have a solid solution with Ti ions occupying substi- 
tution positions. 

3. Experimental Details 

The spec tromet er • used for the measurement v^/as a 
conventional delayed coincidence set up c ompr ising of two 
2" X 2" Nal(Tl) scintillators coupled to RCA-56AVP photo- 
multiplier tubas, kept at an angle 135°. The fwhm of the 
prompt curve was 6 nsec w it h the gates set at ~70 keV energy 
in both the channels.. The source to detector distance was 
kept at" 7 cms. The block diagram _of the electronics (mostly 
an ORTEC system) is given in Fig. (2). 
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A magnetic field of 7 kOe was applied to polarize 
the host matrix, perpendicular to the plane containing the 
detectors. A small electromagnet was used for this purpose 
(Atomic Laboratories Inc. U.S.A.). The magnetic field was 
measured accurately with a Hall probe, which in turn was 
calibrated using proton resonance. The photomultiplier 
tubes were adequately shielded from stray magnetic fields 
using a series of Netic-CoHetic magnetic shields [ 4]. 


4. Measurement 


Approximately 16 k counts were collected for each 

44 

of the field directions. A typical time spectrum for 5cCo 
in an external polarizing field is given in Fig. (3). The 
Fourier transforms of the reduced counting rate R(t), 
defined in Chapter II, are plotted in Fig. (4). 


The measured hyperfine field H is given by the 

me as 


relation 


H 


meas 


->true -»• -V 

H , + H, „ - DM^ + H 

ext hf s L 


as defined in Chapter I, The hyperfine field value after 

correcting for the external and demagnetizing field is 

H, - ScCo (RT) = -58 + 4 kOe. The field is not corrected for 
ht — ~ 

Lorentz field [ 5 1. 


The sign of the hyperfine field was determined by 
the usual method, i.e. from the sense of rotation of the 
angular correlation pattern. The approximate uncertainties 
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are, 4% in the frequency determination from the experimsntal 
data, 2fo in the magnetic moment value and in the time 
calibration, which gives an error of 7% in the measured 
hyperfine field. 

5. Discussion 

After the present measurement was published, Colley 

et al , [6] reported hyperfine field measurements on ScCo 

using TDP AC technique. They however have reported a series 

of values for the hyperfine field, ranging from -91 kOe to 

+13 kOe, depending upon the method they used for source 

preparation. They prepared their samples by melting cobalt 

44 

plates containing evaporated Ti activity, covered by thin 
layers of Li and Ca. Li and Ca were used as reducing agents 
to prevent any internal oxidation of Ti to TiD^ [7 1. The 
positive field was attributed to an internal oxidation of 
Ti-^-TiO^. This was further confirmed, when they annealed 
their sample 5, which indeed showed a change in sign of 
hyperfine field from a negative to positive. We, however, 
took all the care to prevent any possible internal oxidation 
of the activity and we accordingly did not observe any 
oxidation effects in either of our measurement on ScF_e and 
ScC_o. Benski and Hohenhemser [S] prepared their samples in 
hydrogen atmosphere and in vacuum for the hyperfine field 
measurement on Sc in Ni and they did not observe any oxida- 
tion effects either. 
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The observed amplitude is ~ 2.5/4. The observed large 
width (6 r^) of the absolute Fourier transform may be from 
1 ) HCP struoture of the cobalt host with unresolved 
quadrupole interactions, 2) all the 5c ions are not in the 
same well defined lattice positions resulting in unresolved 
hyperfino frequency components. However, further experiments 
are necessary to understand this problem. 

Should we expect a smooth variation of the hyperfine 
field values around 5c impurity in the cobalt host , the 
expected hyperfine field for 5cC^ from the systematic s [•9'] 
is = -60k0e, which compares well with the present 

experimental value of -58 kOe. 

The dilute impurity hyperfine fields in the ferromag- 
netic hosts are often assumed to be proportional to their 
host moments at least when there are no loc aid-moments . 

The hyperfine field on 5c (3d 4s ) in Fe and Ni hosts have 
been reported to be -129+4,k0e [ 2] and -25.7+.0.3 kOe [8] 
respectively. The ratios of the reduced hyperfine fields 

(H j,/U, ) in C o and Ni' hosts, normalized to the reduced 
ht h 

hyperfine field in Fe host, turns out to be 0. 8:1.0 
indicating thereby that the fields are proportional to the 
host moments to within 20 %. Later in Chapter V, it has 
been shown that such a conjecture, as expected, is valid 
only very approximately and indeed the hyperfine field 
depends on many other contributions arising due to impurity 
host interactions and due to impurity itself. 
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1 2 

Scandium has a configuration 3d 4s . The orbital 
angular momentum for a transition impurity like 5c in a 
cobalt host is. expected to be largely quenched and there- 
fore the orbital contribution to the hyperfine field may 

be neglected to a first approximation. Because of the 

,12 

configuration, 3d 4s , the d— moment contributions are 
expected to be small. Therefore the major contribution 
to the hyperfine field should come from the conduction 
electron polarization. Sc with its 3d 4s configuration 
would have rather a small contribution due to core 
polarization, because the magnitude of local d-moment in 
the beginning of the transition series are expected to be 
s m al 1 . 

With the recently available accurate data, we have 
recalculated the polarization 'p’ defined in Shirley et al's 
model [ 5 ] (discussed in detail in Chapter V) and the values 
obtained are 0.026 for Fe host, 0.02 for Co host and 0.023 
for Ni host. Using these values and the free atomic 
hyperfine field value of 1 MOe, we get the conduction 
electron polarization contribution = - -p p = -35 kOe.’ 

Therefore, the contributions from the core polarization is 
~ -23 kOe. Using Freeman and Watsons [ 1 0'] estimates for 
the hyperfine field per polarized d-spin (H^^/3d spin = 

-125 kOe), local d-moment at the scandium impurity site 
turns out to be Q-2^g 
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The atomic volume of Sc is 18 K and accordingly 
from the Stearns model, one would expect a large positive 
contribution to the hyperfine field from the volume overlap 
effect, which is not observed. Therefore this model is not 
successful in explaining the experimental field [11]* 

Finally Balabanov and Delyagin [ 1 2] , approximated the 
observed hyperf ine field across a period to a parabola of 
the t yp e , 

[-2.43 + 0.113 (V- 9)^] 

For scandium, V= 3 and Z =18. This model predicts 

0 

a value +113 kOe, much too different both in magnitude and 
sign of the observed field. 

The model due to Shirley et al. therefore seems to 
explain the experimental value both in sign and magnitude. 
However, a further check of this model would be to measure 
the local d-moment using neutron diffraction techniques. 

The volume overlap effect suggested by Stearns and the 
model due to Balabanov and Delyagin are not successful in 
explaining the observed field. It is, however, quite clear 
that conduction electron polarization is a major contributor 
to the field on Sc in Co host. 
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B. HYPERFINE F lELD ' ME A5URE MENT5 ON TANTALUM IN COBALT 
AND NICKEL HOSTS 

1 . I ntroductian 


The spin 
Ta w as used 
and Co hosts, 
in the B "dec ay 
decay scheme is 


rotation of 10.8 ns, 5/2*^, 

to measure the hyperfine f 

The 133-482 keV cascade in 
181 

of Hf was used for this 
given in Fig. (5) [13j. 


482 keV state of 
ield of Ta in Ni 
^^"^Ta, populated 
purpose. The 


2. Source Preparation 

1 8 1 

The 42.4 day, Hf activity was produced through 
an ( n ,Y ) reaction of Hf metal powder (99.9%). The neutron 
irradiation was carried out at the CIRUS reactor, ' B ARC (India). 
F or g factor measurements, active Hf was dissolved in 26 N ■ 
hydrofluoric acid. In this form the quadrupole interactions 
are expected to be small. The liquid' sample containing HfF^ 
was then sealed in a thin walled polythene tube and was u^d 
as such for the measurements. 

The T aCo sample was prepared by a radio frequency 
induction melting of cobalt plates (99.99%) with traces of 
finely powdered Hf activity in an inert argon atmosphere. 

The molt was gradually cooled to room temperature. The 
resulting shining ball was used as such for the measurements.. 

The Tafy. sample was prepared by the rmally diffusing 
the ^®Va activity into a thin nickel foil (99.99%). The 








radioactive Hf in the form of HfCl^ (prepared by heating 
the metal Hf activity in a chlorine atmosphere) was sand- 
wiched between Ni foil and then sealed in a quartz tube in 
vacuum ( 10 cm) . The diffusion was carried out at 
IIQO'^C for 00 h. The active foil was then thoroughly 
cleaned with distilled water to ensure a complete removal 
of the undissolved surface activity and was then used as 
such for the measurements. Another sample of T aN i prepared 
by induction melting was also studied to see the source 
preparation effects. The procedure adopted for the sample 
preparation was similar to the one used for T aCo ^ 

3. Experimental Details 

The spectrometer consisted of two 2” x 2" Nal(Tl) 
scintillators coupled to 56 AVP photomultipliers. The 
detectors were kept at an angle of 135°. The system had 
a fwhm of -“2.5 nsec for the energy gates set at 133-482 keV 
energies, T';;, g f^'C’car measurement was carried out in an 
external polarizing field of 7 kOe. The hyperfine field 
measurements on T a Ni was also carried out using the same 
set up in an external polarizing field of 2 kOe. The 
other experimental details are discussed in part A. 

The spectrometer for the hyperfine field measurements 
on cobalt host consisted of 1/2" x 2" NaI{Tl) and 1" x 1" 
lead loaded plastic (Pilot B) scintillators. The system 
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had a fwhtn resolution of ~ 1 .7 nsed for the 133 — 482 keV 
energy gate settings. The measurements was carried out 
using randomly oriented domain technique (Chapter II, 
eqn. (17)), with the two detectors kept at an angle 180°. 

4 , Mb asu re men t 

Typical TDPAC time spectra obtained for g-factor 
measurements in an external field of 7 kOe is given in 
Fig. (6) and the reduced counting rate R(t) is given in 
Fig. (7). The TDPAC time spectra obtained for T aN i samples 
(both induction melting and diffusion samples) along with 
the corresponding reduced counting rates R(t) and C(t) 
respectively are given in Figs. (8) and (9). Finally, the 
TDPAC spectra along with the reduced counting rate C(t) 
obtained for T aCo sample is given in Fig. (10). Figure (11) 
gives the absolute Fourier transform plot of C(t) for T aCo . 
s amp le . 

The g-factor of 482 keV state is measured to be 
g = 1. 3+0.1. The hyperfine field value in the nickel host 
after correcting for the external and the demagnetizing 
fields is T aN i (RT, Diffusion )= -77+5 kOe and Ta 

Ni (RT, Induction melting) = -75_+ 5 kOe. For TaC_o, the 
data reduction was done by forming the ratio C(t), defined 
earlier. As shown earlier in Chapter II, the TDPAC spectra 
in a randomly orientdi domain sample is in general modulated 
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by two frequencies and 2Wj^. In the present study, 

however, the 2w^ modulation could not be observed because 
the corresponding time period of modulation (72^ ) was 
comparable to the fwhm resolution of the system. From 
the Fourier transform, which indicated two distinct and 
broad fre q uency d is tribut ions , the hyperfine field value of 
T aC o (RT) is deduced to be = (-) 347+.7 kOe. The sign 
is not measured in this experiment, but is expected to 
be negative, based on systematic trends. The basis of error 
estimation has been the same as detailed in part A. 


5. Discussion 


In the g-f actor measurement only one Larmor oreces- 
sion frequency could be observed. Our value is in agree- 
ment with the previously reported value of g = 1.29 ± 0.02 
by Matthias et al , [141 and by Agarwal et al . [15], The 

observed amplitude is 12%, which is close to the value 
expected from the unperturbed correlations. 

Hyperfine field on Ta in Ni has been very extensively 
studied [16-23], The results reported so far along with 
the present values are detailed in Table I. A restudy was 
however undertaken with the objective of probing into the 
composite nature of the field in the Ni matrix, as was 
pointed out by Bozek et al.[22]. Further Vanderledeen [17] 
reported an amplitude of 0.3^, while Cameron et al . [18] 

obtained an amplitude 1,5%. The measurements reported by 
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Bozek et al . [ 22], indicated a steep fall in the amplitude 
of R(t), which lead them to conclude that the hyperfine 
field has a composite nature, much in contrast to the 
observation by Cameron et al . More recently Livi et al . [23] 
reported their moasu rerne nt on T aN i using an induction 
melting sample. 

We in the present measurement studied two sanples 
made by diffusion of Hf activity and by induction melting. 

The sample prepared by diffusion did not indicate any 
observable attenuation in R(t). The method followed by 
us was exactly similar to the one employed by Bozek et al , 
except that we carried out the diffusion for 80 h and 
they did it for 10 h only. Although we did not observe 
any fall in the amplitude of C(t), the Fourier transform 
did indicate a broad frequency distribution with its 
centroid around 150 MHz. It therefore appears that longer 
diffusion times are necessary to get an appreciable frac- 
tion of Hf atoms in substitutional lattice sites. The 
measurement with T aNi induction melted sample also indicated 
a similar behaviour. It, therefore, based on the present 
observations, is concluded that all the observed discrepancies 
may be explained in terms of rather low solubility of Hf 
ions in the Ni matrix, which in turn results in a variety of 
lattice sites occupied by the Hf impurity and thus results 
in a broad spectrum of hyperfine fields observed by sources 
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prepared by different methods. The implanted sample used 
by Cameron et al . with very dilute Hf concentration' did 
indicate a unique hyperfine field, showing that most of 
the Hf ions came to rest in substitutional lattice sites 
after implantation. This may lead us to conclude that 

we do not yet completely understand the source preparation 
problems and' each system should be treated independently. 

Concurrent to the present measurement in cobalt 
matrix, Soares et al.[24 ] have recently reported the 
field on T aC o to be T aC o (RT) = 362.4 (5.0), employing 

e^-y TDPAC. Their sample was prepared by isotope separator 
implantation. It is however interesting to note that in 
both the measurements two distinct frequencies were 
observed, which as has been pointed out by Soares et al. 
allows for the possibilities of two different lattice 
sites acquired by Hf ions in cobalt host. The almost one 
to one correspondence in the results obtained by contrasting 
methods of source preparation indeed supports such a point 
of view. From the intensities of the different frequency 
components it can be concluded that ~30-4Q/ra of Hf ions do 
acquire unique lattice sites. The weak sattellite peaks 
in Fourier transform show that the Hf ions occupy different 
lattice sites as well. However, sharper frequency distri- 
bution observed by Sores et al corroborate well with the 
results of Cameron et al . indicating that that in implantaf- 
tion, most of the ions come to. rest at unique lattice sites 
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(probably substitutional). It, therefore, appears that 
for Hf in Co and l\l i hosts ion implantation method is more 
favourable compared to metallurgical methods. 


Expecting a smooth variation in the hyperfine field 
values around Ta impurity, the expected field on Ta, from 
the systematics (Fig. (1-3); Chapter V) are -325 kOe and 
-75 kOe for the cobalt and nickel hosts respectively. 

These values compare reasonably well with the experimental 
values. A comparison of the ratios of the reduced hyperfine 
fields (H|_^_p/Vi^) in the cobalt and nickel hosts normalized 
to the reduced field in Fe host indicates that the fields 
are proportional to the host moment only within 30/o. 

The orbital contribution for a transition series 
.32 

impurity like Ta ( 5d 6s ) is expected to be largely quenched. 
The major contribution for such a transition series is 
expected to arise mostly through conduction electron 
polarization (Hj-^p). The other significant contributions 
are the core polarization field (Hj,p) and the volume 
overlap polarization (H^^), 


Using Shirley et al . ’ s prescription [ 5], the conduction 

electron polarization is expected to be Hj-^p (Ta Co ) = -290 kOe 

and Hj,g.p ( TaNi) = -110 kOe, where the free electron hyperfine 

field for Ta is taken to be H = 8.3 MOe . Such a H__p, 

ns CEP 

contribution leaves behind a contribution of -60 kOe 
(cobalt host) and +30 kOe (nickel host) to bo explained by 
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a sum of core polarization and volume overlap polarization. 

The atomic volume of Ta is 10.3 ^ and therefore a finite 

volume overlap contribution may be expected. Using our 

recent estimates presented later in Chapter V, the volume 

C o 

overlap contribution turns out to be = +72 kOe and 

IM i 

= +53 kOe. This in turn would imply a core polarization 

r*' IV t * 

field of kOe and Hj^p = ~25 kOe, for the two 

hosts respectively. Using Freeman and Watson's [10 ] 
estimates of H^^p = 1200 k0e/5d, unpaired spin, the deduced 

core polarization fields indicate that Ta impurity in 
cobalt and nickel host develop a d-moment of ~ 0,1 y and 
~0.02Ug in cobalt and nickel hosts respectively. However, 
if we do not take volume overlap contribution the expected 
d-moments are obviously smaller. Such small values of 
d-moments are indeed supported by the temperature dependence 
of the hyperfine field measurements on T aNi , reported by 
Barret et al. [ 25]. More sensitive and direct probes like 
diffuse neutron scattering experiments vjhic h directly yield 
the d-moment value and hence H|-p are expected to be useful 
in determining the validity of various hyperfine field 
models and as well as in understanding the detailed breakdown 
of various contributions to the hyperfine field on Ta in 
CoandNitoosts. 

The Balabanov and Delyagin [12] model predicts 

T aC o = -295 kOe and H TaNi = -96 kOe in good agreement 

ht — 

with the measured value. The model however does not offer 
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any microscopic details on the origin of the hyperfine 
fields. 

From the above measurement, it could be concluded 

t h at 

(1) Conduction electron polarization gives the dominant 
contribution to the hyperfine field on Ta in Ni and Co, 

(2) All the observed discrepancies on Ta l^i could 

be understood on the basis of the low solubility of 
Hf atoms in l\ii host. 

A similar behaviour is also seen from the measurements 


on Ta in Co host. 
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chapte/L IV 


HYPERflNE INTEPACT70M STUPIES ON 
nmUlTES ANP UAGNETITE OPES 
USING MOSSBAUER EFFECT 


Here, we discuss the Mossbauer effect studies on 
natural iron ores viz. ilmenites (FeTiO^) and magnetites 
(Fe^O^). Mossbauer effect studies on mine ralo gical samples 
is in te res tin g , mainly for two reasons; (1) presence of suf- 
ficient amount of iron in most of iron ores of the earth 
and hence quantitative estimates of iron could be obtained 
using Mossbauer effect, (2) the ability of Mossbauer effect 

to distinguish between the different oxidation states of 
n ■ r- 2+ 3 + 

Fe, i.e., e.g., Fe and Fe > which in turn reflects on 
the geochemical genetics and the subsequent weathering 
history of the sample under perview. 


* Published in Phys. 5tat. Solidi ( a) 321 (1 974). 
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A. ILMENITE5 

1 . Introduction 

Ilmenite is common in nature as an accessory mineral 

in basic igneous rocks, notably the gabbross. Its structure 

is similar to that of hematite ( 86282^5 except that the 

octahedral voids are occupied by Ti^"*" and Fe^"*" ions [ll. 

It is also an ubiquitous mineral in beach sands, concentrated 

by weathering process. The beach sands of Kerala (India) 

and Florida (U,5.A.) are important members of this category. 

During weathering and concentration, the alteration process 

includes oxidation and leaching. The resulting structure is 

disordered and resembles pseudo-brookite [ 2 ], however, the 

arrangement of oxygen atoms as well as their relative 

composition is preserved in this process. Ilmenites from 

3 

beach sands have dominant Fe compared to others [3l. Such 

3 *4* 

an enrichment of Fe ^ indicates the extent of weathering 
the mineral has undergone. From a comparison of measured 
transmission through the sample, with the transmission of 
another sample of known stoichiometry, it is possible to 
make an estimate of iron content of the sample. 

2. Experimental 

The details of the Mossbauer spectrometer and its 
performance are given elsewhere [4]. A Nuclear-Data 512 
multichannel analyzer in a multiscaling mode was used for 
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data acquisition. A Xe-CH^, proportional counter with an 
fwhm of 14% was used for the detection of 14.4 keV -rays. 

The samples used V'jere finely powdered and the absorbers 
were prepared by sandwiching the samples between thin paper 
sheets. All the spectra were recorded at room temperature. 
Since we used a sinusoidal drive, the spectra were linearized 
by fitting the data to a sine function with the help of a 
computer program. 

3. Measurement 

The Mossbauer spectra obtained for ilmenite samples 
from Kishengarh, Bihar, Kodaikanal, Kerala and IMellore 
(India) and from Florida (USA) along with a spectrum of a 
synthetic sample (Temp-Press Inc., USA) are given in 
Figs. (1) and (2). The Mgssbauer parameters are detailed 
in Table I. The errors are the uncertainties in locating 
the peak positions. The statistical uncertainties are 
somewhat less. 

Quantitative estimates of iron was obtained using 
the areas under the peaks. The recoil-free fractions were 
assumed to be same in all the cases. The samples were also 
chemically analyzed and the results are also presented in 
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4. Discussion 

Various phases present in the samples under study 
may be identified using the fact that the different oxida- 
tional states of Fe have different isomer shift. The 
isomer shift was shown to be 

5 = ^ [ I ^.(0) 1^ - I 1,(0)|^n < 

The isomer shift is proportional to the difference 

in the s-electron density at the source and at the absorber 

nuclei. The isomer shift changes by a sizeable amount for 

different oxidation states of Fe. Such a difference arises 

because of difference in the shielding of 3s- electrons by 

3d electrons. Physically, the Ss-electrons spend a part 

of time, farther from the nucleus compared to the 3d-electron. 

The electrostatic potential experienced by them therefore 

depends on the screening effect of inner electrons. Addition 

of a 3d-electron reduces the attractive Coulomb potential 

and this causes the 3s wave function to expand, thereby 

reducing its density at the nucleus. Thus, a removal of 6th 

d-electron in going from ferrous to ferric state, would 

result in an increased charge density at the nucleus. Further 

57 

since 5R/R is negative for Fe, an increase in charge 

density would consequently imply a proportionate decrease 

57p- 

in isomer shift for Fe: i.e. 

(0)I^> I «3+ '*2 + 



88 


Systematic measurements of the Mossbauer spectra of 
various iron compounds coupled to theoretical calculations 
by Walker et al.[ 5 ] and Canon et al , [ 6]. have established 
that the range of values of isomer shifts is distinct for 
various oxidation states of Fe. In such conditions, it 
has been possible to identify the different oxidation states 
of iron from the measured isomer shift and this has been the 
basis of present investigation. Table II below gives 
typical ranges of the isomer shifts referred to sodium 
nitropursside standard [7 1. 


Table II 



^ 2 + 

3 + 

4 + 

6+ 

Oxidation state 

Fe 

Fe 

Fe 

Fe 

Isomer shift 

+ 1 . 0- 

+0.4- 

+0.3 

- 0.6 

(mm/sec) 

+ 1 . 8 

+0.9 



Typical example 

FeF2 

FePD^ 

SrFeO^ 

K^FeO^ 


The ilmenites chosen for the present study may be 
grouped into two categories, viz., 

i. Ilmenites of primary high temperature origin, and 
ii. Secondary ilmenites of placer origin. 

i. Ilmenites of primary high temperature origin 

Four ilmenites viz., synthetic, Kishengafh, Bihar and 
Kodaikanal fall in this category. The synthetic sample was 
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chosen for its high purity and ideal stoichiometry. The 
Mossbauer spectrum of the synthetic sample is a doublet 
with an isomer shift of +1 . 23- ir, m/sac (w.r.t. SNP), which is 
typical of Fe^"*". The quadrupole splitting is 0.71 mra/sec. 
These values are in fair agreement with the earlier reported 
values [8]. The existence of two small satellite peaks 
reported by Ruby and Shirans [ 8 ] is also confirmed in the 
present studies. Those peaks may be due to atmospheric 
oxidation of ferrous iron to a slight extent. The isomer 
shift of the ilmenite samples from primary igneous associa- 
tion, viz., Bihar, Kishengarh and Kodaikanal, cluster around 
1.3 rp.m/sec, showing that these samples have a near stoichio- 
metric composition, with iron mostly in Fe state. The 
observed quadrupole splittings are also nearly the same as 

that observed for the synthetic sample. The possible 

3 “4^ 

oxidation of these samples leading to Fe is seemingly 
low and an upper limit of is obtained for Fe^"*" out of 
the total Fe content. 

ii. Secondary ilmenites of placer origin 

The major source of these ilmenites are beach sands 
(placer) deposits, which have been formed by accumulation 
of heavy minerals (with higher specific gravity) due to 
continuous weathering and wave action. The samples from 
the beaches of Kerala, Florida and Saidapuram fall into 
this category. The Mossbauer spectra of Kerala and Florida 
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and Saidapuram fall into this category. The Mossbauer 

spectra of Kerala and Florida samples clearly indicate the 

coexistence of Fe'~'*' and Fe ions. The dominance of Fe 

in both these spectra clearly establish that these samples 

have undergone considerable vi/eathering and that the Florida 

sample has undergone relatively through more cycles of 

2 + 3 + 

weathering. The ratio of the Fe to Fe for ilmenite ore 
of Florida origin was measured previously [2] using the 
Mossbauer effect and our result is in agreement with theirs. 

In addition the Saidapuram sample shows a doublet 
superimposed on a magnetic hyperfine structure pattern, which 
is typical of a solid solution (1-x) FeTiO^-Cx) Fe202’ 
spectra thus revealed that the sample contains a substantial 
fraction ( 9Q^o) of hematite with about 10% ilmenite. 

Thus the ilrnenites of primary high temperature origin 
could be unambiguously distinguished from the secondary type 
on the basis of isomer shift. Further, the relative 
intensities of Fe"' and Fe in the sample indicates the 
geochemical genetics and as well as the subsequent weathering 
the sample has undergone. The quantitative estimates 
obtained using the Mossbauer spectroscopy are given in 
Table I. From these data one finds that these values are 
comparable to the values obtained from the chemical methods. 
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B. MAGNETITES 
1 . Introduction 

Magnetite is one of the most abundent oxide minerals 
of iron and has its occurrence in igneous and in sedimentary 
sources. In igneous rocks it is commonly of late stage 
magmatic seggregation origin and is associated with tha 
basic magma. Because of its high temperature genesis, 
titanium enrichment in its composition is frequently 
observed. For magnetites in sedimentary formations, many 
theories have been put forth [9-13], 

Since the conditions of formation of magnetites 
differ, it is expected that compositional variations should 
also exist. As pointed out earlier, a study of such 
variations, would enable us to understand the genetic 
history of the sample. Structurally, magnetite is an 
inverse spinel. From its structure, composition for this 
mineral should contain 68.97% Fe^O^ end 31,03% FeO. However, 
this ratio varies in a natural sample, which may be inter- 
preted as due to a variance in oxygen supply during its 
formation. In the present studies two samples from Bihar 
(India) (one of igneous origin and the other from banded 
iron formations), and two samples (one from Karnataka and 
another from Tamilnadu (India)), associated with banded 
iron formations [14] were taken up. For comparison, one 
synthetic sample (Temp-Press Research Inc. USA) was also 


stud ied . 
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2. Experimental 

The experimental conditions were essentially similar 
to those described in part A. 

3 . Measure ment 

The Mossbauer spectra obtained for these magnetite 
samples are given in Figs. (3) and (4). The Mossbauer 
parameters are provided in Table III. The errors are the 
errors in the peak positions. The statistical uncertainties 
are however somewhat less. 

4. Discussion 

The magnetites chosen for the present study may be 
classified into two categories, viz., 

i. Magnetites of high temperature origin 
ii. Magnetites of possible low temperature origin 

i. Magnetites of high temperature origin 
a. Synthetic magnetite 

The spectrum given in Fig. (3(a)) is clearly a super- 
position of two iron spectra with different isomer shifts 
and different hyporfine fields corresponding to sites A 
(tetrahedral) and sites B (octahedral), having ferric and 
ferric + ferrous ions respectively. The analysis was done 
by separating the spectra following the procedure suggested 
by Bauminger et al, [ 15], The present values compare well 
with the earlier reported values [15-17] . 
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Table II. Mossbauer Parameters of Magnetites of Different 
Geochemical Origins. 


S.No. 

Sample 

Mossbauer Parameters 


S ite 

Isomer 

shift* 

( mm/ sec ) 

Hyperf ine 
f iel d 
( kOe) 

Quadrup ole 
sp litt.in gs 
( mm/ sec ) 

1 

Synthetic Magnetite 

T etra. ( A) 

0,5 + 0. 1 

487+4 

0. 1±0.'i 



Octa. (B) 

1 .1 + 0.1 

441+4 

0. 1 ±0-'* ■ 

2 

Igneous Magnetite 
(Bihar, India) 

Te tra. (A) 

Octa. (B) 

’ 0 . 6 + 0 . 1 

1 .0 + 0.1 

483 + 4 

453+4 

0. 1 + 0.1 

0 • OHh 0 • 1 

3 

Magnetite, BMQ 
(Karnataka, India) 

Te tra. ( A) 

Octa. (B ) 

0.5+0. 2 

0. B±0.2’ 

4 97 + 8 

480+8 

0 . 0+0 . 2 

0 . 1 + 0 . 2 

4 

Magnetite, BMQ 
(Tamilnadu, India) 

Te tra. ( A) 

Octa. ( B ) 

0 . 5j^0 . 2 

0.9+0. 2 

497+8 

476+8 

0.1+0. 2 

0. 1+0.2 

5 

Magnetite, BHQ 
(Bihar, India) 

Te tra. ( A) 

Octa. (B) 

0.4 +0.2 

0. B+0. 2 

526+8 

508+8 

0 . 0+0 . 2 

0 . 2+0 . 2 


* w.r., t. Sodium Nitro Pr us side. 



b. 


Magnetites from igneous associations of Bihar 
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The similarity of Mossbauer parameters with those 
obtained for the synthetic sample indicate a near stoichio- 
metric composition of the sample, which is expected from a 
sample of high temperature origin. Further, there is a 
clear indication of a quadrupole doublet. Fig. (3(b)), which 
is possibly due to the titanium enrichment in the sample. 
This sample was chosen from an igneous association in 
nature and microscopic investigations in fact indicated the 
presence of thin lamallae of titanium-rich portions along 
the cleavage planes in magnetite. These observations 
further confirmed the titanium enrichment, as evidenced in 
the Mossbauer spectra. 

ii. Magnetites of possible low temperature origin 

The samples from Karnataka, Tamilnadu and Bihar fall 
into this category. Of these, the samples from Karnataka 
and Tamilnadu are from banded magnetite silica formations, 
and the Bihar sample is from banded haematite quartzite 
formations. The Tamilnadu sample is, however, an enrichment 
of magnetite within these banded magnetite formations due to 
a near complete leaching of silica. 

The spectra of these three samples - Fig. (4) are 
once again a superposition of two spectra. The isomer shift 
obtained, however, indicated the presence of iron in ferric 
state only. Further the magnetic hyperfine field on the 
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two sites were found to be almost equal. These observa- 
tions suggest a possible oxidation of magnetite to hematite 
These banded magnetite quartzite formations are typical 
metamorphosed horizons and although initially magnetite 
was one of the dominant iron oxide mineral, much of it 
got oxidized to hematite during metamorphism, which thus 


explains the 

indication of Fe' 

3+ . 

xn 

the 

spectra. A com 

pa 

r ative 

study 

of these 

spectra 

with 

that 

of 

the synthet 

ic 

s amp le 

clearly brings 

out the 

near 

abse 

nee 

of ferrous 

ir 


phase, thus supporting above explanation. 

The sample taken from the banded hematite-quartzite 
formation of Bihar also indicated the presence of a ferric 
phase, which again should be due to near complete oxidation 
of magnetite. 

Thus it is possible to differentiate magnetites of 
igneous origin from those of low-temperature origin. In 
the former a near stoichiometric structure is evidenced, 
while in the latter, a dominance of ferric phase over 
ferrous phase occurs. Indication of titanium enrichment 
in the magnetites of high-temperature origin could be very 
clearly observed in the spectrum. 
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Vi LUTE. IMPURITY HYPERPINE PIELVS IN 
FERROMAGNETIC Fe Co ANV NX MATRICES 
-THEIR ORIGINS ANV COMPARISONS OF 
THE BEST EXPERIMENTAL VALUES 
WITH THE EXISTING MOVELS 


A. IIMTRDDUCTIDN 

Samoilov at al . [1], first reported that non-magnetic 

impurities dissolved in ferromagnetic hosts like Fe, Co and 
Ni experience large magnetic hyperfine fields. Since then, 
the magnitude and the sign of these dilute impurity hyper- 
fine fields have been measured on almost all the impurities 

for the iron host. These fields have been very useful in 

\ 

understanding a variety of solid state and nuclear physics 
phenomena. The hyperfine fields have proved to be very 
useful, as a diagnostic device in providing a local probe 
of the electronic state of an atom in metals. These large 
internal hyperfine fields, which are not otherwise available 
in the laboratory are extensively used in the study of the 
magnetic properties of very short-lived nuclear excited 
states [2] . . ' 
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Precise theoretical estimates of these fields are 
practically impossible,, due to the complexity of the 
Hamiltonian and the many body techniques involved. The 
fact that these fields show a very regular systematic 
variation, when plotted as a function of impurity elec- 
tronic configuration, prompted semi-theoretical approaches, 
to decode the basic mechanisms involved. The present 
status seems to be that we understand all the major 
mechanisms responsible for these fields [3-14], however a 
precise breakdown of their relative strengths is yet to 
be achieved. 

The present chapter describes our studies on the 
systematica of these fields. Comparisons of the experi- 
mentally measured hyperfine fields, with the predictions 
of the various models are presented. 

B. HOST MOMENT DEPENDENCE OF THE HYPERF.INE FIELDS 

For a detailed study of the host moment dependence 
and for the model comparisons, all the experimental data 
reported so far on the dilute impurity hyperfine fields in 
Fe, Co and Ni hosts was collected and collated [I5l. From 
this tabulation, the 'best' experimental values were chosen, 
after a careful consideration of various factors like, the 
impurity concentration, the temperature of the measurement 
and the experimental technique enployed. These values are 
tabulated in reference [161 . A plot of these values vs th_ 
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impurity electronic configuration is given in Figures (1-3). 

The magnitudes of the hyperfine fields for an 

impurity in Fe, Co and Ni hosts are often assumed to be 

proportional to the magnetic moments of the respective host 

matrices = 2.22 = 1.72 and = 0.604 Pg) . 

However, this assumption was investigated in greater detail. 

To examine this, the 'reduc ed fields* (H,_/u.) for Fe, Co 

hf h 

and Ni hosts were calculated. A plot of these ratios, i.e. 

impurity electronic configuration is given in Fig.i'4). It 
is clearly seen that the observed hyperfine fields are not 
strictly proportional to the host moments. The deviations 
are usually larger than 10^. 

Such deviations are in fact expected, should we 
consider a detailed dependence of the hyperfine fields on 
the host and impurity moments. Thus, based on Kushida et al . ’ s 
discussion [17 I, we have, 


hf ±mp 


3 d-4s 


core 

3i 


• ^h 


Where p ^ and pi are the impurity and the host moments 
respectively. a represents the contributions due to core- 

polarization of the impurity shells by the impurity moment, 
ns ■ 

a is the 'self* polarization of the conduction electrons, 
due to impurity moment, a^*^ represents the contributions 
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c o x* o 

due to an admixture of 3d-4s wave functions. a^ 

o s 

represents core polarization due to host and is the 

conduction electron polarization due to the host moment. 

All these contributions are discussed in a somewhat greater 
detail in the following section. It can be very well 
realized that it is rather futile to expect any rigid 
proportionality with the host moment as has been very 
clearly evidenced in Fig. (4). 


C. HYPERFINE FIELD MODELS : COMPARISON WITH EXPERIMENTAL 
VALUES 

1 . Hamiltonian 


A general Hamiltonian for an impurity atom, charac- 
terized by the quantum numbers L, S, and J, embedded in a 
ferromagnetic host can be given as [4], 

H + 2lig H^^.S +xt.5 -Yj H^-l” + 2vgYi < >{( -k) *5 ^ 


* 2 V'l " ^ * 2y Uj 

X r . 

1 



E 

i 


^ ^ 

3 ( r . *3 . ) ( r . • I ) 
1 11 



( 2 ) 


where is the Bohr magneton, Y -r = nuclear 

D i i N 

g-f actor and U., is the nuclear magneton. V . is the crystal 

IM C t 

field potential, describing the influence of surrounding 

atoms. H is the magnetic exchange field due to host 
ex 

moments. XL*S represents the electronic spin-orbit 



interaction term. 
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represents the contributions from 

the spin polarization of the conduction electrons via a 

Fermi contact interaction [1 8 1 . kS-I represents the core 

spin-polarization of the filled shells of the impurity 

atoms and this also generates a Fermi contact hyperfine 

field. k is a parameter measuring the extent of core 

polarization, l.*I represents the hyperfine interaction 

between the magnetic field set up due to the orbital motion 

"fc h 

of the electrons, 1. being the angular momenta of the i 
electron with a position vector r^. The last term represents 
nuclear d ip ol ar -el ec tron spin dipolar interaction term, the 
external, demagn etizing and the Lorentz fields have been 
discussed in Chapter I. 

The first three terms in the Hamiltonian are largest 
and thus determine the zeroth order wave function. The 

-I** ^ 

relative magnitudes of and X L*5 allow us to distingurah 

between two cases and accordingly a further discussion could 
be carried over under two subheadings: 

-»■ 

2. Strong Crystal Fields - Won rare-earth Solutes 

3. Weak Crystal Fields - Rare-earth solutes (XL-S'3>\/ ) 

C T 




2. Strong Crystal Fields - Non rare-earth Solutes {XL*S>>V ) 

C, T 

F pr large crystal fields, the orbital angular momentum 
is almost quenched, i.e., <L^>« 0. This is largely true f.or 
transition series solutes. For example, in iron series the 
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unquenched orbital angular momentum contribution is 
typically ~ 10 kOe, which is only a few percent of the 
free atomic hyperfine field. Such small residual values 
of angular momenta, arise essentially due to the distortion 
of the orbital wave function by the weak spin-orbit part. 
The Hamiltonian is further simplified as the dipolar term 
also vanishes, if the ferromagnetic host has a cubic 
symmetry (Dipolar terms vanish for 5-state ions as well). 
Under such conditions, the Hamiltonian reduces to 

" - ''af * (3) 

Under these conditions, the crystal fields determine the 
zeroth order electron wave function. The exchange term 
can then be [4] used to get the spin expectation value of 
the ion in the lattice. The last two terms can be combined 
to define an effective hyperfine field at the nuclear site, 
i. e , , 

Yj + 2ygYj <r“^ >{k) Sa (4) 

i.e. + 2Ug <r"^> (k) "s (5) 

Hence the effective magnetic hyperfine field in non rare- 
earth solutes can then be thought to arise largely from 



Conduction electron polarization 

: H|.g.p 

ii. 

Core polarization 

: Hj,p, and 

iii. 

Overlap polarization 



and 
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i. Conduction electron polarization : 

In this case, one considers that the magnetic centres 
act as spin dependent scattering potentials, that strongly 
perturb the spatial distribution and the energy of the free 
conduction electrons described by 

^l<r = 4!(a) (6) 

The spin-dependent scattering potentials perturb the 
up and down spins differently and produce an oscillatory 
spin polarization. Such spin polarized conduction electrons, 
then generate hyperfine field at an impurity atom at a 
distance R via a Fermi contact term. 


8 , 2 2 2 

H|-^p = — Iu^(R)I >[iU;_j^{0)| -1 i|»^(Q)j 3 (7) 

R R 

Such an oscillatory spin polarization may arise from (a) 
(a) Friedal's charge perturbation [8-9], {b).the RKKY type 
moment perturbation [19l and (c) the s-d interband mix- 
ing [20] . 

The conduction electron polarization could be given 
as [1 2 ] , 


’^CEP “ 

where Hj represents 
electron polarization 
polarization) and H 

s 


H_ + H (8) 

S s 

a summed contribution to the conduction 
by all the neighbours (transferred 
is the polarization induced by the 
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impurity atom (self-polarization) due to its local d{f) 

■moment. 

a. Daniel and Freidal model and Campbell extension 

An impurity in a matrix is considered as a spherically 
symmetric potential V, the depth of which is adjusted to 
yield the right screening charge, which in turn is determined 
by Friedal's sum rules and the phase-shifts. In the first 
order, this potential depth is determined by the charge 
difference AZ, between the impurity and the host 
(after properly taking care of the closed shells and 
d-scree ning ) . For iron 1, [21], while for impurity it 

will be the number of valence electrons. The width of the 
potential is decided by the number of conduction electrons 
per atom. 

The magnetic exchange coupling between the host 
moment and the conduction electrons spin splits the conduction 

* 

A small correction to this valance difference was proposed 
by Blatt [ 22 3 , who suggested that the introduction of a 
impurity results in an outward displacement of the host 
neighbour atoms, which in turn is equivalent to a removal 
of some fraction of ionic charge from the region occupied 
by the impurity, i.e. 

SV *5 y 

A = AZ - where ~ is the fractional change 

in the size of the cell occupied by the impurity. 
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band resulting in a difference of energy for the two spins 
(up and down). Thus the conduction electrons with dif.ferent 
spins (up or down) see different scattering potentials i.e, 

= V - e; = V + e and = 2e (9) 

where e is the s-d exchange energy. This results in an 
unequal scattering for the up and, down spins. Such an 
unequal scattering results in a difference in the radial 
wave function of the two spin types ( rfi and ) leading 

to a net spin polarization. The conduction electron hyper- 
fine field is then given by the proportionality, 

E 2 2 

^CEP “ / I - 1^4, (0) I n^(E)]dE (10) 

o * 

Daniel and Friedal assumed, one s-type conduction 
electron per host-atom. A ferromagnetic coupling was 
assumed to explain the observed experimental data correctly. 
The model was quite successful because it could correctly 
predict the observed change in sign of the hyperfine fields 
in ' sp’ solutes for AZ = 3 (Sb-5n). 

Campbell [10] later improved the Daniel and Friedal' s 
model. He showed that ~ 0.25 5-type conduction electrons/host 
atom is a more realistic number. This as shown in Fig, (5) 
resulted in very different estimates of CEP compared to 
Daniel and Friedal model. Campbell assumed an antiferro- 
magnetic s-d coupling, whereas Daniel and Friedal model 
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Conductian electron .polarization at the iRpurity 
site predicted by ,a )' Daniel and Friedal '•s'-Hodei' 
b) C u3!fp u;*.! ij. *a fiiodel'^ plotti-^d as of-unction of 
volums .Z. of ,i( 5 purity otouu ■ ' 
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implies a ferromagnetic coupling. Campbell also calculated 
the self-polarization contribution for the transition 
impurities by modifying the impurity potential from, 

V,- v^ = -2e 

to a difference, 

V . - y , = -2e + Jy . (11) 

T r 1 

where J appropriately modifies the potential to take care 
of impurity moment. Campbell got the conduction electron 
contribution as 

H^g-p = A(z)P^ = A(z) (0.6 + 0.4 (12) 

U and being the impurity and host moments respectively. 
This model predicts a sharp decrease in the hyperfine field 
values with an increase in the charge difference. Further 
a largo AZ results in V » e(the s-d exchange energy) and 
accordingly one would expect y and V_|^to be not too different 
and hence the resulting and are also similar. In these 
cases the field is largely decided by the mean spin polariza- 
tion of the conduction band. On the other hand, when 
V - e i.e. A Z ~ 0, as is generally true for transition 
series solutes, the difference in spin up and spin down 
wave functions is significant and the resulting hyperfine 
field is therefore critically dependent on the s-d exchange 
energy and is opposite to the mean spin polarization. 



Such a model though intituitively simple, however, 
suffers from a large handicap in providing exact ’numbers’ 
especially due to a complicated band-structure and Fermi 
surfaces in alloys, and one has to rely heavily on some 
model dependent parameters and band structure calculations. 

The other polarization mechanism is the RKKY [ 19 I 
type interaction, where the oscillatory spin density of 
the conduction electron arises due to a coulomb exchange 
interaction between the localized d-electrons and the 
conduction electrons. Such an interaction always results 
in a net positive conduction electron polarization. 

Finally, s-d interband mixing aris,es due to the fact that 
in real transition metals, there is no clearcut distinction 
between the s- and d-elecirons [20] and one therefore 
considers an effective exchange interaction J' through 
interband mixing of the conduction and the local moment 
electron orbits. Such an interaction in most cases leads 
to a negative CEP. 

b. Model due to Shirley et al . 

Another successful approach to estimate conduction 
electron polarization hyperfine field was initiated by 
Shirley and Westenbarger [3 ] and later improved upon by 
Shirley, Rosenblum and Matthias 14]. Without going into 
the detailed dynamics of the origin of these fields, these 
authors proposed that would be a scaled fraction of th 



free atom hyperfine field i.e. the hyperfine field generated 
by an unpaired outermost s-electron. The model however 
explicitly assumed a host moment proportionality. Thus, 

= p • P. • H (13) 

CEP ^ h ns 

where the factor p-absorbs all the host and solute dependent 
factors like say, the electron transfer in alloying and the 
difference in s-electron density at the nucleus in a free 
atom and in the metallic lattice. Inspite of the simplicity 
of the model, it has provided a reasonably reliable and an 
extremely simple vjay to calculate the conduction electron 
polarization. 

The value of p is obtained from the experimental 

1 0 

data for dilute impurities with a filled core, i.e. nd 
structure (Cu, Zn, Ag, Cd , Au, Hg) , where the conduction 
electron contribution is expected to be predominant. Using 
more recent data, the values of p are recalculated and are 
tabulated in Table I, It may however be worthwhile to note 
that this model does not give any z-dependence of p, in 
contrast to D an iel-Friedal-C ampbell ' s model, where the 
conduction electron contribution was found to be critically 
dependent on the electronic configuration of the - impurity . 

c. Stearns experiments: Scaling rules 

Stearns [12] was able to measure the conduction 
electron polarization due to impurity neighbour shells (Hj.) 
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Solute 

H 

ns 

(MOe) 

Host 

Fe 

Host 

C o 

Host 

Ni 

H^^(kOe) 

P 

H^^(kOe) 

p 

H^^(kOe) 

P 

3d'’ °Cu 

2.7 

-21 3 

0.078 

-1 57 

0.058 

- 45 

0.016 

4d"’ °Ag 

5.0 

-447.2 

0.089 

-108 

0.021 

-122.4 

0.024 

4d'' °Cd 

6.5 

-340 

0.052 

-268 

0.041 

— 6 8.6 

0.013 

5d''°Au 

19 . 8 

-1170 

' 0.059 

-797 

0.042 

-272 

0.013 

5d''°Hg 

25.0 

-440 

0.017 

-370 

0.014 

- 80 

0.003 

Average 

’ P ' 


0.059 


d.035 


0.013 
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0.027 


0.021 
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in the Fe host| VJith a careful analysis of the Mossbauer 
and the pulsed NMR results on ordered Fe^Si alloy, she was 
able to experimentally deduce the total conduction electron 
polarization of all the neighbour shells. Assuming the 
additivity of contributions of different shells, she gave, 


Fe 


I M. &H. 


[\!=1 


where gives the number 

is the hyperfine field 
is h 

the i shell. Assuming 
host moments, one obtains, 
These values could also be 
hyperfine field values of 
expected to be dominant. 


= -145 kOe (14) 

is h 

of iron atoms in the i shell, 

contribution from an Fe atom in 

that is proportional to the 

H = -120 kOe and = -46 kOe 

I i 

obtained independently through 
Cu, Au and Ag where CEP is 


ii. Core polarization 

The notion of core polarization was first invoked by 
Sternheimer [23], who used it to explain the observed large 
negative hyperfine fields on transition series solutes. The 
core electrons are paired and henge do not contribute 
to the hyperfine field in the first order. In the second 
order, the exchange interaction between the Sd-electrons 
of the solute and the solute core electrons results in an 
unequal spatial distortion of the radial wave function for 
the up and down spins, leaving a net spin density at the 
nuclear site. The usual convention is that the magnetic 
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hyperfine field is positive, when it is parallel to the 
direction of magnetization. Since positive field corres- 
ponds to a positive spin, the spins are said to be positive 
when antiparallel to the direction of external field. The 
3d spin of the impurity orients opposite to the direction 
of the external magnetization. This in turn 'attracts* 
the core ns + (n = 1 ,2,3) electrons and repel nsi electrons, 
leaving thereby a net negative spin density at the nuclear 
site and hence a negative field. Figure (6) brings out 
this contention explicitly. Freeman and Vi/atson [7] have 
employed unrestricted Hartree-Fock method to calculate the 
contribution to the hyperfine field from each of the core 
ns shells. Their calculations show that out of the core 
ns shells, the 2s shell's contribution is maximum. Freeman 
and Watson also calculated the core polarization hyperfine 
field/unpaired d-electron to be as -125 k0e/3d spin, 

-375 k0e/4d unpaired spin and -1200 k0e/5d unpaired spin. 
These fields/spin are expected to be constant throughout 
a transition series, which has been verified experimentally. 
Thus, in Shirley's model, the total field is given as, 



V.pH +2<S>H 


ns 


nd 


u + H 

"cep CP 


(15) 


Thus a localized d-moment at the solute atom which in turn 
polarizes its own core ns shells, can give large negative 
core polarization fields. The neutron diffraction data [24] 
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'by the 3d, electro'fjs, the ‘arrows ind ic ate. 'the,- average , , 
■'.positions of ''the 'inner, ou'tex'. and'.:3c.i ' ei-8ct.ron.„spinS' '..and, , 

', ■' . the, effects excHa.nge , xnte.rac'tipn., : ■.■',' '■,. ■.: 
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in fact has given adequate evidence for the existence of 
such local d-moments . Accordingly, the d-moments derived 
from equation (16), using the experimental hyperfine field 
data- may be compared with the measured d-moment values from 
neutron diffraction experiments. For such a comparison, 
from the values of p given in Table I, the CEP estimates 
were calculated. The CEP estimates thus obtained were 
subtracted from the experimental hyperfine fields to get 
the core polarization contribution to the hyperfine field. 

From these estimates of the core polarization contributions, 
and using the calculations of Freeman and Watson [6,7] the 
localized d-moment values were calculated, and are tabulated 
in Table II and. plotted in Fig. (7), along with the experi- 
mental values obtained from neutron diffraction measure- 
ments, wherever applicable. 

Recently Cheng-jyi Song et al. [25] have reported a 
first direct measurement of the 2s-3s spin densities using 
Mossbauer effect with electron spectroscopy. Their experi- 
ment, to some extent, verify the calculations of Freeman 
and Watson [5-7], Duff and Das [26] and of Liu and Bagus [27] . 
The satisfactory predictions of the CP model in giving 
localized d-moments, as confirmed by diffuse neutron 
scattering data, along with Cheng-jyi Song’s experimental 
confirmation of theoretical calculations do give us a 
reasonable amount pf confidence to conclusively say that 
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T able 1 1 


Series 

Element 


Host Fe 

Host Co 

Host 

Ni ■ 

C ale . 

•lig Expt.yg 

Calc.yg Exp t. -Jig 

Calc .pig 

Expt . lig 

3d 

Sc 

0.4 


0.2 

0.1 



Ti 







■ V 

0.2 


0.1 




Mn 

1 .0 

0 + 0.5 

0.6 

2.4 



Fe 

1 . 8 

2.22 

2.1 

1.9 



Co 

1.3 

2. 1+0.3 

1 .2 

0.8 



Ni 

0.7 

1 .2 + 0.3 

0.2 

0.7 



Cu 

0,4 


0.5 

0.1 


4d 

Nb 

0.3 


0.3 

0.0 



Mo 

0.2 


0.2 




T c 

0.4 


0.4 

0.0 



Ru 

0.8 

0, 7+0.4 

0.8 

0.5 



Rh 

0.9 

0 « 5^0 . 3 

0.7 

0.4 


• 

Pd 

0.9 

0 . 3 + 0 . 3 

0.3 

0.4 



Ag 

0.4 



0.1 


5d 

Ta 

0.2 


0.1 

0.0 



W 

0.1 


0.1 




Re 

0.1 


0.0 




0 s 

0.3 

0.1+0. 3 

0.4 

0.1 



Ir 

0.4' 

0 . 2+0 . 3 

0.4 

0.2 



Pt 

0.2 

. 0.1+0. 3 

0.2 

0.1 



■ Au 



0.1 

0.1 
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the CP is a dominant contribution in the transition series 
solutes. Such a model further strengthens Shirley's model 
and his predictions vis a vis CEP estimates. 

iii. Overlap polarization 

The other approaches in calculating transferred 
hyperfine fields is the chemical bonding approach [13]. 

In this the impurity atom is treated as a ligand coordinated 
to its nearest magnetic neighbours. As a result of overlap 
and/or covalency, the ligand's outermost electrons get 
spin polarized, subsequently generating a hyperfine field 
via a Fermi contact interaction. Following such an approach 
Stearns [12], suggested the possibility of a volume overlap 
term and Shirley [13] tried to explain the fields especially 
on d'^'^s^p'^ solutes by calculating the overlap of 3d host 
and 5s impurity wavef unctions . The contact field resulting 
due to an overlap polarization is essentially positive [7] . 

a. Stearns Volume Overlap Model 

This model postulates a misfit volume overlap contri- 
bution due to valence ns-like electrons, which remain bound 
to the impurity site. Such an overlap polarization, therefore, 
occurs only if the volume of the , impurity exceeds the effec- 
tive volume available to it in the host matrix. Thus, when 

V. > V, the valence ns-like electrons, overlap with the 

imp host’ 





host matrix and get positively pgj,^j,zed by an ajpount 
proportional to the volume misfit of the volume of the 
solute atom, Steasn? gave a simple lineajr relation fop 
such a polarisation, i,e., 

p^{z) = C{V^ - V^) (16) 


where C is a constant characteristic of the hogt, is 

the volume of the solute and V is the effective volume 

o 

available to the impunity on removing one host atpjn. The 


predictions of this mqdei were calculated both by using 

equation (17) ,and by using the experimental data as ris. 5 cribed 

P 0 

below. Using Stepprits estimates of = -145 kOe , 

Co {\l i 

Hj = -117 kOe and -44 kOe and the scaling rules 

given by her, i.e. 


(17) 


the conduction electron polarization contribution wqs 
calculated. These values were then subtracted from the 
experimental hyperfine field values of the solutes, which 
do not develop a moment (i.e., non transition and non pare 
earth series solutes) to get volume overlap contribution. 

It is thus implicitly assumed ihat the total contribution 
to hyperfine fislds in such solutes ( d s p type) arisps 
predominantly from the condi^otion electron polarization 
and the volume overlap polarization. The volume part 
so obtained was then divided by the appropriate free atomic 
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hyperfine field to get the experimental volume overlap 

polarization values These are plotted as a function 

of the impurity volume in Fig. (S). A least squares 

fit of these values to the equation (16) gave both the 

effective volume V and the constant C, for the three hosts 

o 

as shown in Fig. (s). These least square fitted values were 
then used to calculate the predictions of the model using, 
the relation (16). The impurity volumes were taken from 
Ref. [12,28]. The calculated volume overlap polarization 
p^ values are plotted along with the experimental values 
as a function of impurity Z, in figure (9) and are tabulated 
in Table III. The general gross features for the non- 
transition series solutes may be reproduced by this model. 

As expected, the overlap of the solute s-electrons with 
the polarized d-electrons of the host can give significant 
contributions to the hyperfine field. This contribution 
may be to a crude approximation given by the volume overlap 
effect. Although the gross features may be explained by. 
this model, as one would expect there are large deviations 
in some cases and these deviations may probably be explained 
by invoking the microscopic structure of the ion. 

b. Shirley's Overlap Model 

Outlining the limitations and the inadequacies of 

the Daniel and Frieda!' s model for large values of solute 
1 Q 3 n 

valency (d s p , A Z ») , Shirley [ 1 3 ] suggested that the 
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large fields in these solutes arise due to 5s shell polariza- 
tion by the host 3d-spins. Such a polarization could come 
either via a P auli-dist ortion effect and/or by a covalency 
term. Orthogon alization of the host 3d spins and the 
impurity 5s spins lead to an admixture of the Id-character 
into 5s-orbital with spin parallel to the 3d spin. This 
results in a Fermi contact hyperfine field s.t. 

;• = !< 3d I 5s>l^ (18) 

The amplitude <3dl5s> givesthe extent of overlap. Such an 
overlap of atomic wave functions is shown to explain the 
observed trends of the hyperfine fields in 5sp solutes. 

iv. Balabanov and Delygin's model 

Balabanov and Delyagin [14] studied the dependence of 
the magnetic hyperfine fields on the number of electrons in 
the outer shells of the impurity atom (i.e., o.n the total 
number of (n-1)d electrons, ns electro-ns and np electrons 
where n = 4 , 5 and 6 respectively for the elements of the 
period IV, V, VI). They approximated the observed regularity 
by means of a parabola, that is symmetrical about the middle 
of the period, i.e., 

= a + b (‘v - 9)^ (19) 

where V is the number of outer electrons and a and b are 
constants. The model assumes that the electrons in the outer 
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shell polarizes the inner filled shell to give rise to a 
negative hyperfine field on the impurity. The positive 
contribution predominates whenjv - sj ^5. The magnitude 
of the positive contribution is also determined by the 
polarization of the internal shells, but does not depend 
on the value of the effective magnetic moment of the outer 
electrons of the impurity. Such a polarization is caused 
by a direct interaction of the electrons in the filled 
shells and the polarized electrons of the conduction band 
of the matrix. 

Obviously enough, the predictions of this model 
depends on how best the experimental data could be fitted 

to a parabola. Also the assumption that the maximum 

1 .2 

negative field is obtained for minimum value of jv - 9 \ , 

i.e., in the middle of the period is not tenable, as the 
maximum is obtained for three different impurities for the 
three hosts (Figs. 1-3). ' 

This model unfortunately does not offer any micro- 
scopic explanation of the polarization except that the 
outer shell electrons contribute to the hyperfine field. 
The experimental data presented in Figs. (1-3), show that 
the data can only be approximately represented by means of 
a parabola. In Fig. (10), the least squares fitted curve 
obtained from this model (fitted to equation (19)) along 
with the experimental data is given. The predictions of 
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this model were compared with the experimental values for 
the fourth group in Fe, Co and Ni hosts and in 20 cases, 
the deviation is more than 20%. In some cases, the sign 
of the field predicted by this model is also not in agree- 
ment with the reported experimental sign. The host moment 
proportionality assumed by the model, has already been seen 
to have a limited validity. Further, many of the systematic 
trends of the observed hyperfine fields are not explained 
by this model [1 6]. 


3. Weak Crystal Fields-Rare Earth Solutes (XL*S >>y _) 

cf 


For weak crystal fields the Hamiltonian reduces to 
the form, 


-»• -► 


-3 
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Symbolic ally , 


H, . ^ H . . . . + H , + H , + H + H 

hr spxn-orbit exch orb core cond 


+ H 


spin-dipolar 


( 20 ) 


Further discussion could be carried out for Lanthenide 
(4f) and Actinide (5f) solutes separately. 
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i. Lanthenide solLrtes. 

a. Orbital contribution 

The 4f solutes retain much of theix fyee ion character, 
due to a very effective shielding of the 4f electrons by tins 
5s and 5p shells and accordingly the magnitudes of the 
observed field is primarily decided by the large unquenched 
orbital angular momenta of the 4f shell. The contribution 
due to the 4f electron orbital momenta is practically an 
order of magnitude larger than other contributions. The 
orbital contribution is expressed as, 

= 2Ug < 4f { r"^ I 4f>< j{f N|j J >J (21) 

The reduced matrix elements I1 n}} are tabulated by 

-3 

Elliot and Stevenson [29] and the <4f | r~ j 4f > values 
have been determined both [30] by theory and experiment, 

b. Core polarization contribution 

The mechanism of core polarization in the rare earths 
is the same as in the transition series solutes, discugsed 
earlier except for the fact that here the 4f electrons are 
responsible for core polarization. Assuming core polarization 
to be constant, within the 4f seiies , Freeman and Watsqn [ 7] 
gave an empirical relation based on available data, that the 
core contribution could be given as 

H 


c ore 


9Q (Qj - 1 ) J (kOe) 


( 21 ) 
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However with a closer scrutiny and more accurate 
calculations by the same authors indicated that the field/ 
unpaired spin was not constant and did vary considerably 
within the 4f series. 

c. Conduction electron polarization 

An estimate of conduction electron polarization could 
be obtained by Shirley's model. This contribution is 
however expected to be small owing to smaller 6s densities 
at the nucleus. 

Thus the fields on rare earth (4f) solutes is primarily 
decided by the orbital part, which more or less keep their 
free atomic like character. The other contributions are in 
most cases an order of magnitude smaller. The sign of the 
field (especially a change in the sign near the middle of 
the 4f shell) can be qualitatively understood assuming an 
antiferromagnetic coupling between the 3d electron spins 
and the 4f spin, which in turn, align antiparallel (parallel) 
to the angular momenta in the first (second) half of the 4f 
shell[6b3. 


ii. Actinide solutes 

In actinides [31] as compared to the lanthenides. the 
5f electrons are not that well shielded as the 4f electrons 

and in these elements, the spin orbit and crystal fields 
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are of comparable magnitudes resulting in partial quenching 
of the angular momenta. Further due to a significant overlap 
of 5f wave function with 6d, 6s, Is wave function and 
accordingly configuration mixing effects become important 
for such solutes. Under such conditions, .the core polariza- 
tion effects become more important. Dunlop [32] has given 
on the basis of experimental data, an empirical relation 
for core contribution as, 

H = -(560 + 100) (g,-1) J (kOe) (22) 

core “ J 

The experimental hyperfine field me asurejnents for 
actinides have been quite limited. The reasons being 
experimental difficulties arising due to! 1) a-decay after 
effects, 2) existence of more than one charge states, and 
3) the daughter nuclei are also radioactive.. 

With limited data available, no studies on systema- 
tics could be undertaken. 

D. CONCLUSIONS 

From the plots . of the best experimental dilute 
impurity hyperfine fields for Fe, Co and Ni matrices, many 
interesting systematic trends could be seen. Comparison of 
the ’reduced hyperfine field' values for Fe, Co aqd Ni 
hosts showed that the usual assumption, that the impurity 
hyperfine field values are proportional to the host moments. 
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is found to be not strictly valid and the deviations 
usually larger than 10%. Because of the ^lomplexity of the 
Hamiltonian a|id the many body techniques involved, the 
detailed calculatiops to obtain the dilute impurity hype^ffine 
fields in f er:|omagnetio hosts do not eeeinto be feasible i,n 
the near future. Therefor®, most of the existing models 
for the dilute impurity hyperfine fields are semiempir ical 
and are derived from meagre theoretical estimates and from 
the observed syst,amatics of the fields. From the present 
comparisons of the experimental hyperfine field values with 
the predictipns of the various existing models, it is 
concluded that the core-polarization aqd the conduction 
electron polarization model suggested by Shirley et al. is 
reasonably satisfactory for the transition series solutes. 

The volume overlap model due to Stearns can give the general 
gross features and crude estimates for the nsp solutes. The 
model due to Balabanov and Delyagin also gives the general 
tre nd s . 

Thus, the dilute impurity hyperfine fields is one area,, 
where the experimeni,s are far ahead of the theoretical calculs 
tions. More systematic and rigorous microscopic calculatipps 
for the different possible contributions are necessary to 
further understand the mechanisms involved regarding the 
origin of the dilute in^urity hyperfine fields. 
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